THE SPECTRAL CURVE OF THE EYNARD-ORANTIN RECURSION VIA 

THE LAPLACE TRANSFORM 
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Abstract. The Eynard-Orantin recursion formula provides an effective tool for certain enumer- 
ation problems in geometry. The formula requires a spectral curve and the recursion kernel. We 
present a uniform construction of the spectral curve and the recursion kernel from the unstable 
geometries of the original counting problem. We examine this construction using four concrete ex- 
amples: Grothendieck's dessins d'enfants (or higher-genus analogue of the Catalan numbers), the 
intersection numbers of tautological cotangent classes on the moduli stack of stable pointed curves, 
single Hurwitz numbers, and the stationary Gromov-Witten invariants of the complex projective 
line. 
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1. Introduction 

What is the mirror dual object of the Catalan numbers? We wish to make sense of 
this question in the present paper. Since the homological mirror symmetry is a categorial 
equivalence, it does not require the existence of underlying spaces to which the categories 
are associated. By identifying the Catalan numbers with a counting problem similar to 
Gromov-Witten theory, we come up with an equation 

(1.1) x = z + - 



as their mirror dual. It is not a coincidence that ( 1.1 ) is the Landau-Ginzburg model in one 
variable [2] 00]. Once the mirror dual object is identified, we can calculate the higher-genus 
analogue of the Catalan numbers using the Eynard-Orantin topological recursion formula. 
This recursion therefore provides a mechanism of calculating the higher-order quantum 
corrections term by term. 

The purpose of this paper is to present a systematic construction of genus spectral curves 
of the Eynard-Orantin recursion formula |25|, I27| . Suppose we have a symplectic space X 
on the A-model side. If the Gromov-Witten theory of X is controlled by an integrable 
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system, then the homological mirror dual of X is expected to be a family of spectral curves 
S. Let us consider the descendant Gromov-Witten invariants of X as a function in integer 
variables. The Laplace transform of these functions are symmetric meromorphic functions 
defined on the products of S. We expect that they satisfy the Eynard-Orantin topological 
recursion on the B-model side defined on the curve S. 

More specifically, we construct the spectral curve using the Laplace transform of the 
descendant Gromov-Witten type invariants for the unstable geometries (g, n) = (0, 1) and 
(0, 2). We give four concrete examples in this paper: 

• The number of dessins d'enfants of Grothendieck, which can be thought of as higher- 
genus analogue of the Catalan numbers. 

• The V'-class intersection numbers (r^ • • • Td„) g n on the moduli space M. g , n of pointed 
stable curves [El [251 EE EE] ■ 

• Single Hurwitz numbers [TJ [231 ES] . 

• The stationary Gromov-Witten invariants of P 1 |62| [66] . 

The spectral curves we construct are listed in Table 1. The Eynard-Orantin recursion 
formula for the single Hurwitz numbers [U [TJ [Ml EZ] and the ^-class intersection numbers 
|25j are known. Norbury and Scott conjecture that the stationary Gromov-Witten invariants 
of P 1 also satisfy the Eynard-Orantin recursion |62j . A similar statement for the number of 
dessins d'enfants does not seem to be known. We give a full proof of this fact in this paper. 



Grothendieck's Dessins 


( 


x = z + \ 
y = -z 


" " " T d n )g : n 




J X = z 2 

\y = -z 


Single Hurwitz Numbers 


( 


x = ze 1 ~ z 
y = e z_1 


Stationary GW Invariants of P 1 


1 


x = 

\y = 


--Z+ 1 - 

z 

-- -log(l + z 2 ) 



Table 1. Examples of spectral curves. 



Let Dg >n (fjLi, . . . , /x n ) denote the weighted count of clean Belyi morphisms of smooth 
connected algebraic curves of genus g with n poles of order (fj,i, . . . , /j, n ). We first prove 

Theorem 1.1. For 2g — 2 + n > and n > 1, the number of clean Belyi morphisms satisfies 
the following equation: 

n 

(1.2) mDg tn (fj,i, . . . , ix n ) = ^(Ml +H -2)D s , n -l(Ml + - 2 > M[n]\{lj}) 

i=2 



Dg-l, n +l(a,P,fJ>[n]\{l}) + D gi,\I\+l( a :t J 'j) D g2,\J\+l(^l i j) 

91+32=9 
IUj={2,...,n} 

where m = (/ij)jg/ for a subset I C [n] = {1, 2, . . . , n}. 
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A),i(2m) = — C m 
2m 

is given by the Catalan number C m = ( 2 ™)- The next case Do,2(Ml> M2) is calculated in 
|44l 145] . Note that the (g, n)-terms appears also on the right-hand side of (1.2). Therefore, 
this is merely an equation, not an effective recursion formula. 
Define the Eynard-Orantin differential form by 

W g D n (t 1 ,...,t n ) = d 1 ---d n Yl Dg, n {^...^ n )e-^ 
where the Wj-coordinates and ij-coordinates are related by 



•+tlnW„) 



tj - 1 tj + 1 



Then 



Theorem 1.2. The Eynard-Orantin differential forms for 2g—2+n > satisfy the following 
topological recursion formula 



(1.3) W g D n (h,...,t n ) 

1 1 

~64 2vri 



i=2 



1 1 \ (t 2 - 1 3 1 , 

+ it 1 r-dh 

t + ti t-hj t 2 dt 



stable 



+ E <W + i(M/)<,|,| + i(-M;) 



91+92=9 

JUj={2,3,...,n} 

This is now a recursion formula, since the topological type (</ , n') of the Belyi morphisms 
appearing on the right-hand side satisfies 

2g' - 2 + ri = (2g - 2 + n) - 1, 

counting the contributions from the disjoint union of the domain curves additively. A 
corollary to the recursion formula is a combinatorial identity between the number of clean 
Belyi morphisms and the number of lattice points on the moduli space M gn that has been 
studied in [lOl [551 [591 E01 E] . 

Corollary 1.3. 



(1.4) r> ff , n ( m ,..., yUn )= e •• e ri 



,>aM=i 



Mi 



Ng tn {2ii - m,--- ,2£ n - fi n ), 



where N gin (/j,i, . . . , fj, n ) is defined by (4- IS). 



The recursion formula (1.3) is a typical example of the Eynard-Orantin recursion we 
discuss in this paper. We establish this theorem by taking the Laplace transform of (1.2). 
This is indeed a general theme. For every known case of the Eynard-Orantin recursion, 
we establish its proof by taking the Laplace transform of the counting formula like (1.2). 
For example, for the cases of single Hurwitz numbers [24^ [57] and open Gromov-Witten 
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invariants of C 3 |8H I82j. the counting formulas similar to (1.2) are called the cut-and-join 
equations [291 EH 113 EE [80] . 

The Laplace transform plays a mysterious role in Gromov-Witten theory. We notice its 
appearance in Kontsevich's work [46] that relates the Euclidean volume of M g , n an d the 
intersection numbers on -M 9jn , and also in the work of Okounkov-Pandharipande [65J that 
relates the single Hurwitz numbers and the enumeration of topological graphs. It has been 
proved that in these two cases the Laplace transform of the quantities in question satisfies 
the Eynard-Orantin recursion [10, 24, 26, 55, 57J for a particular choice of the spectral curve. 

Then what is the role of the Laplace transform here? The answer we propose in this 
paper is that the Laplace transform defines the spectral curve. Since the spectral curve is a 
B-model object, the Laplace transform plays the role of mirror symmetry. 

The Eynard-Orantin recursion formula is an effective tool in certain geometric enumer- 
ation. The formula is originated in random matrix theory as a machinery to compute the 
expectation value of a product of the resolvent of random matrices ([I], [21] )• In |25tl27| Ey- 
nard and Orantin propose a novel point of view, considering the recursion as a mechanism 
of defining meromorphic symmetric differential forms W gn on the product E n of a Riemann 
surface S for every g > and n > 0. They derive in |25|I27| many beautiful properties that 
these quantities satisfy, including modularity and relations to integrable systems. 

The effectiveness of the topological recursion in string theory is immediately noticed 
[23l [5Tj [70] . A remarkable discovery, connecting the recursion formula and geometry, is 
made by Marino [51J and Bouchard, Klemm, Marino and Pasquetti jS]. It is formulated as 
the Remodeling Conjecture. This conjecture covers many aspects of both closed and open 
Gromov-Witten invariants of arbitrary toric Calabi-Yau threefolds. One of their statements 
says the following. Let X be an arbitrary toric Calabi-Yau threefold, and £ its mirror curve. 
Apply the Eynard-Orantin recursion formula to X. Then W gn calculates the open Gromov- 
Witten invariants of X. The validity of the topological recursion of [251I27] is not limited to 
Gromov-Witten invariants. It has been applied to the HOMFLY polynomials of torus knots 
[9], and understanding the role of quantum Riemann surfaces and certain Seiberg-Witten 
invariants [35]. A speculation also suggests its relation to colored Jones polynomials and 
the hyperbolic volume conjecture of knot complements |13j . 

From the very beginning, effectiveness of the Eynard-Orantin recursion in enumerative 
geometry was suggested by physicists. Bouchard and Marino conjecture in [7] that particu- 
lar generating functions of single Hurwitz numbers satisfy the Eynard-Orantin topological 
recursion. They have come up to this conjecture as the limiting case of the remodeling con- 
jecture for C 3 when the framing parameter tends to oo. The spectral curve for this scenario 
is the Lambert curve x = ye~ y . The Bouchard-Marino conjecture is solved in [H 1241 [57] . 
The work [24] also influenced the solutions to the remodeling conjecture for C 3 itself. The 
statement on the open Gromov-Witten invariants was proved in [11} l8T| 182] , and the closed 
case was proved in (5J [H3] . 

The Eynard-Orantin topological recursion starts with a spectral curve E. Thus it is 
reasonable to propose the recursion formalism whenever there is a natural curve in the 
problem we study. Such curves may include the mirror curve of a toric Calabi-Yau threefold 
13 EI] , the zero locus of an A-polynomial [13 ES], the Seiberg-Witten curves [33], the torus 
on which a knot is drawn [S], and the character variety of the fundamental group of a knot 
complement relative to SL(2,C) [13J. Now we ask the opposite question. 



Question 1.4. // an enumerative geometry problem is given, then how do we find the 
spectral curve, with which the Eynard-Orantin formalism may provide a solution? 
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In every work of [HI EH [25l [26l E3 EH EH EQl E21 ED E2] , the spectral curve is 
considered to be given. How do we know that the particular choice of the spectral curve 
is correct? Our proposal provides an answer to this question: the Laplace transform of 
the unstable geometries (g, n) = (0, 1) and (0, 2) determines the spectral curve, and the 
topological recursion formula itself. The key ingredients of the topological recursion are the 
spectral curve and the recursion kernel that is determined by the differential forms Wo,i 
and Wo,2- In the literature starting from [25], the word "Bergman kernel" is used for the 
differential form Wo, 2- But it has indeed nothing to do with the classical Bergman kernel 
in complex analysis. It is also treated as the universally given 2-form depending only on 
the geometry of the spectral curve. We would rather emphasize in this paper that this 
"kernel" is the Laplace transform of the annulus amplitude, which should be determined 
by the counting problem we start with. 

Although it is still vague, our proposal is the following 



Conjecture 1.5 (The Laplace transform conjecture). If the unstable geometries (g,n) = 
(0, 1) and (0, 2) make sense in a counting problem on the A-model side, then the Laplace 
transform of the solution to these cases determines the spectral curve and the recursion 
kernel of the Eynard-Orantin formalism, which is a B-model theory. Thus the Laplace 
transform plays a role of mirror symmetry. The recursion then determines the solution to 
the original counting problem for all (g, n) . 



The Eynard-Orantin recursion is a process of quantization [9. 35j. Thus the implication of 
the conjecture is that quantum invariants are uniquely determined by the disk and annulus 
amplitudes. For example, single Hurwitz numbers hg jfl are all determined by the first two 
cases ^o,(/Lti) an d ^0,(^1,^2)- The present paper and our previous work [231 E7| establish this 
fact. The Lambert curve is the mirror dual of the number of trees. 

The organization of this paper is the following. In Section [2] we present the Eynard- 
Orantin recursion formalism for the case of genus spectral curve. Higher genus situ- 
ations will be discussed elsewhere. Sections [3] and [4] deal with the counting problem of 
Grothendieck's dessins d'enfants. We present our new results on this problem, which are 
Theorem and Theorem 1.2 We are inspired by Kodama's beautiful talk [Hj (that is 
based on |45j) to come up with the generating function of the Catalan numbers as the spec- 
tral curve for this problem. We are grateful to G. Gliner for drawing our attention to |44j . 
The counting problem of the lattice points on M. g , n of [HU ESI EHJ [60] is closely related to 
the counting of dessins, which is also treated in Section [4j The Eynard-Orantin recursion 
becomes identical to the Virasoro constraint condition for the f/j-class intersection numbers 
on M. g ,n- We discuss this relation in Section [5J using Kontsevich's idea that the intersection 
numbers on A4 g>n are essentially the same as Euclidean volume of M.g >n - Section [6] is de- 
voted to single Hurwitz numbers. In our earlier work |24| I57j we used the Lambert curve as 
given. Here we reexamine the Hurwitz counting problem and derive the Lambert curve from 
the unstable geometries. We then consider the Norbury-Scott conjecture [62J in Section u\ 
which states that the generating functions of stationary Gromov-Witten invariants of IE" 
satisfy the Eynard-Orantin recursion. We are unable to prove this conjecture. What we 
establish in this section is why the spectral curve of [62j is the right choice for this problem. 

The subject of this paper is closely related to random matrix theory. Since the matrix 
model side of the story has been extensively discussed by the original authors [27], we do 
not deal with that aspect in the current paper. 
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2. The Eynard-Orantin differential forms and the topological recursion 

We use the following mathematical definition for the topological recursion of Eynard- 
Orantin for a genus spectral curve. The differences between our definition and the original 
formulation found in \25\ [2~T] are of the philosophical nature. Indeed, the original formula 
and ours produce the exact same answer in all examples we examine in this paper. 

Definition 2.1. We start with P 1 with a preferred coordinate t. Let S C P 1 be a finite 
collection of points and compact real curves such that the complement E = P 1 \ S is 
connected. The spectral curve of genus is the data (£,vr) consisting of a Riemann surface 
£ and a simply ramified holomorphic map 

(2.1) 7r:S9ti — > vr(t) = x £ P 1 

so that its differential dx has only simple zeros. Let us denote by R = {p±, . . . ,p r } C £ the 
ramification points, and by 

U = U r j=1 Uj 

the disjoint union of small neighborhood Uj around each pj such that it : Uj — > Tr(Uj) C P 1 
is a double-sheeted covering ramified only at pj. We denote by i = s(t) the local Galois 
conjugate of t £ Uj. The canonical sheaf of £ is denoted by fC. Because of our choice 
of the preferred coordinate t, we have a preferred basis dt for K, and d/dt for /C _1 . The 
meromorphic differential forms W g ^ n {t\, . . . , t n ), g = 0, l,2,...,n = 1,2,3,..., are said 
to satisfy the Eynard-Orantin topological recursion if the following conditions are 
satisfied: 

(1) WV(t)etf (£,/C). 

(2) W , 2 (h,t 2 ) = (f^ffs - 7r*(^gfs G #°(E x E,/C® 2 (2A)), where A is the diagonal 
of E x E. 

(3) The recursion kernel Kjfah) £ H°(Uj x C, (/C^ 1 (8) /C)(A)) for t £ Uj and t x £ C is 
defined by 

(2 - 2) Kj{tM) - 2 KMH^l 1 

The kernel is an algebraic operator that multiplies dti while contracts (it. 

(4) The general differential forms W g>n (ti, . . . , t n ) £ H° (E n , K,(*R)® n ) are meromorphic 
symmetric differential forms with poles only at the ramification points R for 2g — 
2 + n > 0, and are given by the recursion formula 



(2.3) W g<n (t l ,t 2 ,...,t n ) = ^-.J2 ( f Kj(t,h 



Wg-l >n +l(t, t, t 2 , ■ ■ ■ , t n ) 



No (0, 1) terns 

+ J! W giM + 1 (t,t I )Wg 2!l J l + 1 (i,tj) 

91+92=9 
IUj={ 2 ,3,...,n} 

Here the integration is taken with respect to t £ Uj along a positively oriented simple 
closed loop around pj, and tj = (tj)i e / for a subset I C {1, 2, ... , n}. 
(5) The differential form W\^{ti) requires a separate treatment because ^0,2(^1,^2) is 
regular at the ramification points but has poles elsewhere. 



(2.4) Wl , 1 {t 1 ) = ^Y.l Kj(t,ti) 



ttt / \ * dx(u) ■ dxiv) 

W 0t2 (u,v)+TT* 



(x(u) — x(v )) 5 



u=t 
v=t 
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Let y : £ — >-C be a holomorphic function denned by the equation 

(2.5) W ,i(t) = y(t)dx(t). 

Equivalently, we can define the function by contraction y = ixWo t \, where X is the 
vector field on £ dual to dx(t) with respect to the coordinate t. Then we have an 
embedding 

£9 1 1 — > (x(t),y(t)) £ C 2 . 



Remark 2.2. The recursion (2.3) also applies to (g,n) = (0,3), which gives Wo,3 in terms 
of Wo,2- In [25, Theorem 4.1] an equivalent but often more useful formula for Wq^ is given: 

W ,2(t,t 1 )Wo,2(t,t2)W 0t2 (t,t 3 ) 



(2.6) Wb,3(tl,t2,t3) = I 

™ 3=1 JU 



dx(t) ■ dy{t) 



3. Counting Grothendieck's dessins d'enfants 

The A-model side of the problem we consider in this section is the counting problem of 
Grothendieck's dessins d'enfants (see for example, |71| I72j) for a fixed topological type of 
Belyi morphisms [3]. Gromov-Witten theory of an algebraic variety X is an intersection 
theory of naturally defined divisors on the moduli stack A4 gtn (X) of stable morphisms from 
n-pointed algebraic curves of genus g to the target variety X. Since we are considering 
tautological divisors, their 0-dimensional intersection points are also natural. These points 
determine a finite set on M. g - n via the stabilization morphism. If we expect that the 
Gromov-Witten theory of X satisfies the Eynard-Orantin recursion, then we should also 
expect that the counting problem of naturally defined finite sets of points on M g , n may 
satisfy the Eynard-Orantin recursion. 

Pointed curves defined over Q form a dense subset of M g ,n- To specify n, we need to 
use Belyi morphisms. When we identify a curve over Q with a Belyi morphism, a natural 
counting problem arises by considering the profile of the Balyi morphism at the branched 
points. In this way we arrive at canonically defined finite sets of points on .A^^ n . 

More specifically, consider a Belyi morphism 

(3.1) b : C^P 1 

of a smooth algebraic curve C of genus g. This means b is branched only over 0, l,oo £ 
P . By Belyi's Theorem [3], C is defined over Q. Let qi, . . . ,q n be poles of b of orders 
(pi, . . . , fi n ) S . This vector of positive integers is the profile of b at oo. In our enumera- 
tion we label all poles of b. Therefore, an automorphism of a Balyi morphism preserves the 
set of poles point-wise. 

A clean Belyi morphism is a special class of Belyi morphism of even degree that has profile 
(2, 2, ... , 2) over the branch point 1 £ P . We note that a complex algebraic curve is defined 
over Q if and only if it admits a clean Belyi morphism. Let us denote by D g ^ n (m, . . . , [i n ) 
the number of genus g clean Belyi morphisms of profile (/jli, . . . , fi n ) at oo € P . This is the 
number we study in this section. 

We first derive a recursion equation among D g>n (fjLi, . . . , fi n ) for all (g, n). This relation 
does not provide an effective recursion formula, because D s>n (/xi, . . . , fj, n ) appears in the 
equation in a complicated manner. We then compute the Laplace transform 

F g D n ( Wl ,...,w n ) = £Wmi,---,aO e-^+-+^"), 

/xi,...,/j n >0 
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and rewrite the recursion equation in terms of the Laplace transformed functions. We then 
show that the symmetric differential forms 

W g D n = d 1 ---d n F g D n 

satisfy the Eynard-Orantin recursion formula. This time it is an effective recursion formula 
for the generating functions W^ n of the number Dg tn {fM\, . . . , fx n ) of clean Belyi morphisms. 
Grothendieck visualized the clean Belyi morphism by considering the inverse image 

(3.2) r = 6- 1 ([o,i]) 

of the closed interval [0, 1] C P 1 by b (see his "Esquisse d'un programme" reprinted in [72J). 
This is what we call dessin d 'enfant. It is a topological graph drawn on the algebraic curve 
C being considered as a Riemann surface. We call each pre- image of G P 1 by b a vertex 
of r. Since b has profile (2, . . . , 2) over 1 £ IP 1 , a pre-image of 1 is the midpoint of an edge 
of r. The complement C \ T of Y in C is the disjoint union of n disks centered at each qi. 
By abuse of terminology we call each disk a face of V. Then by Euler's formula we have 

2-2 5 = |&- 1 (0)|-|£T 1 (l)|+ra. 

A dessin is a special kind of metric ribbon graph. A ribbon graph of topological type (g, n) 
is the 1-skeleton of a cell-decomposition of a closed oriented topological surface C of genus g 
that decomposes the surface into a disjoint union of 0-cells, 1-cells, and 2-cells. The number 
of 2-cells is n. Alternatively, a ribbon graph can be defined as a graph with a cyclic order 
assigned to the incident half-edges at each vertex. When a positive real number, the length, 
is assigned to each edge of a ribbon graph, we call it a metric ribbon graph. A dessin is 
thus a metric ribbon graph with the same length given to each edge. We usually consider 
this length to be 1, so the distance between and 1 on P 1 is measured as ^. 

The concrete construction of [54J gives a Belyi morphism to any given dessin. Thus the 
enumeration of clean Belyi morphism is equivalent to the enumeration of ribbon graphs, 
where we assign length 1 to every edge. The original interest of dessins lies in the fact that 
the absolute Galois group Gal(Q/Q) acts faithfully on the set of dessins. 

An alternative description of a Belyi morphism is to use the dual graph 

(3.3) f = 6- 1 ([l,ioo]), 
where 

[1, ioo] = {l + iy\0<y<oo}cF 1 

is the vertical half- line on P 1 with real part 1. This time the graph T has n labeled vertices 
of degrees (/Ui, . . . ,/x n ). Since we consider ribbon graphs in the context of canonical cell- 
decomposition of the moduli space A4 Si „, we use the terminology dessin for a graph T dual 
to a ribbon graph T. This distinction is important, because when we count the number of 
ribbon graphs, we consider the automorphism of a graph that preserves each face, while 
the automorphism group of the dual graph, i.e., a dessin, preserves each vertex point- wise, 
but can permute faces. In this dual picture, we define the number of dessins with the 
automorphism factor by 

1 dessin of 
type (g,n) 

where T is a dessin of genus g with n labeled vertices with prescribed degrees (fii, . . . , yu n ), 
and Auto(r) is the automorphism of T preserving each vertex point-wise. 
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Our theme is to find the spectral curve of the theory by looking at the problem for 
unstable curves (g, n) = (0,1) and (0,2). The dessins counted in Z?o,i(m) for an integer 
fx € Z + are spherical graphs that contain only one vertex of degree \i. Since any edge of 
this graph has to start and end with the same vertex, it is a loop, and thus fi is even. So let 
us put fi = 2m. Each graph contributes with the weight l/|Aut£>(r)| in the enumeration of 
the number ZJq.iCm)- This automorphism factor makes counting more difficult. Note that 
the automorphism group of a spherical dessin with a single vertex is a subgroup of Z/(2m)Z 
that preserves the graph. If we place an outgoing arrow to one of the 2m half-edges incident 



to the unique vertex (see Figure 3.1), then we can kill the automorphism altogether. Since 
there are 2m choices of placing such an arrow, the number of arrowed graphs is 2mDo,i (2m). 
This is now an integer. By a simple bijection argument with the number of arrangements 
of m pairs of parentheses, we see that 

1 /2m^ 

(3.5) 2mX> ,i(2m) = C„ 



m + 1 \ m 

where C m is the m-th Catalan number. We note that the Catalan numbers appear in the 
same context of counting graphs in [37] . 



Figure 3.1. An arrowed dessin d'enfant of genus with one vertex. 
Define the Laplace transform of L>o,i(/^) by 



(3.6) = J2 A),i(2mK 



-2mw 

m=l 

Then the Eynard-Orantin differential 



Wg 1 = dFfa = Y1 2mD ,i{2m)e- 2mw dw = - J] C m eT 2mw dw 
m=l m=l 

is a generating function of the Catalan numbers. Actually a better choice is (see 
( . n _f 1 1 1_ , 2_ _5_ 14 42_ 

_ l17 J> l17 ilj tAj iXj iAj 

The radius of convergence of this infinite Laurent series is 2, hence the series converges 
absolutely for |x| > 2. The inverse function of z = z(x) near (x, z) = (oo, 0) is given by 

(3.8) x = z+ 1 . 

z 

This can be easily seen by solving the quadratic equation z 1 — xz + 1 = with respect to 
z, which is equivalent to the quadratic recursion 

Cm+i = ^ ] Cj ■ Cj 

i+j=m 
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of Catalan numbers. To take advantage of these simple formulas, let us define 

(3.9) x = e w 

and allow the m = term in the Eynard-Orantin differential: 

(3.10) ' ' ^ 



w $i = - E c - 



m=0 



^,2m+l ' 

Accordingly the Laplace transform of Do,i(2m) needs to be modified: 

™ oo 1 

-2mw ^ ^ " '~ " ^ 



(3.11) 



w 



2m) 



x 



2 m 



log X. 



m=l m=l 

Although numerically Z?o,i(0) = 0, its infinitesimal behavior is given by 

D 0A (2m) 



lim 

m— 5>0 



r 2m 



which is consistent with 



lim 2mDox(2m) 

m— >0 



- logx 
C = 



1. 



From (3.7) and (3.12), we obtain 
(3.12) 



W D x 



-z(x) dx. 



In light of (2.5), we have identified the spectral curve for the counting problem of dessins 
D 9>n (n). It is given by 

(3.13) ( 



z + 
— z 



To compute the recursion kernel of (2.2), we need to identify ^0,2(^1,^2) for the other 
unstable geometry (g, n) = (0, 2). In the dual graph picture, Do,2(Ml> M2) counts the number 
of spherical dessins V with two vertices of degree fj,± and fi2, counted with the weight of 
l/|Aut£)(r)|. The computation was done by Kodama and Pierce in |45} Theorem 3.1]. We 
also refer to a beautiful lecture by Kodama [44j . 

Proposition 3.1 (|45|). The number of spherical dessinss T with two vertices of degrees j 
and k, counted with the weight of l/|Aut£>(r)|, is given by the following formula. 

(J_ ( 2k \ 

2k \ k) 



(3.14) 



^0,2(^1,^2) 



1 -L. ( 2 A ( 2k \ 

4 j+k \ j )\k ) 



1 ( 2 A( 2k \ 

I, j+k+l \jJ\k) 



tn = 0, [12 = 2k ^ 

= 2j ^0,/j 2 = 2k^0 
/xi = 2j + l,fi 2 = 2k + l 



All other cases L>o 2(^1) M2) = 0. Here the automorphism group Aut£>(f) is the topological 
graph automorphisms that fix each vertex, but may permute faces. 

Remark 3.2. The first case is irregular. For /ix = 0, the second vertex has an even degree, 
and hence we have Ck/(2k) graphs. Note that this graph has k + l faces due to Euler's 
formula 2 = 1 — k + (k + 1). The degree vertex has to be placed in one of these faces, 
which makes the total number of graphs 



k + l 
2k 



1 

2k 



2k 
k 
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However, we are counting only connected graphs. Hence degree vertices are not allowed 
in our counting. 

In general the number of dessins satisfies the following: 



Theorem 3.3. For g > and n > 1 subject to 2g — 2 + n > 0, the number of dessins (3.4) 
satisfies a recursion equation 



(3.15) /xi-D s ,nOl>---,Mn) = ^2(^1 + H ~ 2)D g>n -i (mi + H ~ 2 , M[n]\{l,i}> 

i=2 



+ 



a+^=/ii- 



a(3 



D g -i t n+l{.Oi, f3, M[n]\{l}) + ^si.l/l+lC^M^-^fla^JI+lC^iMj) 



ffi+92=3 
7Uj={2,...,n} 



where \xi = (fJ>i)i£i for a subset I C [n] = {1, 2, . . . , n}. TTte Zasi sum is over all partitions 
of the genus g and the index set {2, 3, . . . , n} into two pieces. 

Remark 3.4. Note that when g\ = and 7 = 0, D g ^ n appears in the right-hand side of 



(3.15). Therefore, this is an equation of the number of dessins, not a recursion formula. 



Proof. Consider the collection of genus g dessins with n vertices labeled by the index set 
[n] = {1, 2, . . . , n} and of degrees (/ii, . . . , fi n ). The left-hand side of (3.15 ) is the number 



of dessins with an outward arrow placed on one of the incident edges at the vertex 1 . The 
equation is based on the removal of this edge. There are two cases. 

Case 1. The arrowed edge connects the vertex 1 and vertex j > 1. We then remove the 



edge and put the two vertices 1 and j together as shown in Figure 3.2 This operation is 
better described as shrinking the arrowed edge to a point. The resulting dessin has one less 
vertices, but the genus is the same as before. The degree of the newly created vertex is 
H\ + — 2, while the degrees of all other vertices are unaffected. 





Figure 3.2. The operation that shrinks the arrowed edge to a point and joins two 
vertices labeled by 1 and j together. 

To make the bijection argument, we need to be able to reconstruct the original dessin from 
the new one. Since both \i\ and \ij are given as the input value, we have to specify which 
edges go to vertex 1 and which go to j when we separate the vertex of degree /ii + fij — 2. 
For this purpose, what we need is a marker on one of the incident edges. We group the 
marked edge and /i, — 2 edges following it according to the cyclic order. The rest of the 
fij — 1 incident edges are also grouped. Then we insert an edge and separate the vertex into 
two vertices, 1 and j, so that the hrst group of edges are incident to vertex 1 and the second 



group is incident to j, honoring their cyclic orders (see Figure 3.2). The contribution from 
this case is therefore 



^(^i + fij - 2)D fl)n _i(/xi + fij - 2, M[n]\{i,j})- 
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Case 2. The arrowed edge forms a loop that is attached to vertex 1. We remove this loop 
from the dessin, and separate the vertex into two vertices. The loop classifies all incident 
half-edges, except for the loop itself, into two groups: the ones that follow the arrowed 
half-edge in the cyclic order but before the incoming end of the loop, and all others (see 



Figure 3.3). Let a be the number of half-edges in the first group, and /3 the rest. Then 
a + /3 = fii — 2, and we have created two vertices of degrees a and (5. 

To recover the original dessin from the new one, we need to mark a half-edge from each 
vertex so that we can put the loop back to the original place. The number of choices of 
these markings is a/3. 






Figure 3.3. The operation that removes a loop, and separates the incident vertex 
into two vertices. 



The operation of the removal of the loop and the separation of the vertex into two vertices 
certainly increases the number of vertices from n to n + 1. This operation also affects the 
genus of the dessin. If the resulting dessin is connected, then g goes down to g — 1. If the 
result is the disjoint union of two dessins of genera gi and 52, then we have g = gi + g%. 
Altogether the contribution from this case is 



a+/3=fi!-2 



Dg-i, n +i(a,P,fJ>[ n ]\{i})+ Dg uW+ i(a,fii)D g2j \j\ + i(P,fJ,j) 



3i+92=ff 
IUj={2,...,n} 



Note that the outward arrow we place defines the two groups of incident half-edges uniquely, 
since one is after and the other before the arrowed half-edge according to the cyclic order. 
Thus we do not need to symmetrize a and (3. Indeed, if the arrow is placed in the other 
end of the loop, then a and f3 are interchanged. 



The right-hand side of the equation (3.15 ) is the sum of the above two contributions. □ 



Remark 3.5. The equation (3.15) is considerably simpler, compared to the recursion for- 
mula for the number of ribbon graphs with integral edge lengths that is proved in |10| 
Theorem 3.3]. The edge removal operation of [10] is the dual operation of the edge shrink- 
ing operations of Case 1 and Case 2 above, and the placement of an arrow corresponds to 
the ciliation of [10]. In the dual picture, the graphs enumerated in [10] are more restrictive 
than arbitrary clean dessins, which makes the equation more complicated. We also note 
that PU Theorem 3.3] is a recursion formula, not just a mere relation like what we have in 
(|3.15|). In this regard, (3.15) is indeed similar to the cut-and-join equation (6.28) of [291177]. 



We will come back to this point in Section [6} 



The relation (3.15) becomes an effective recursion formula after taking the Laplace trans- 
form. 
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4. The Laplace transform of the number of dessins and ribbon graphs 

In this section we derive the Eynard-Orantin recursion formula for the generating func- 
tions of the number of dessins. The key technique is the Laplace transform. 

Since the projection x = z + 1/z of the spectral curve to the x-coordinate plane has 
two ramification points z = ±1, it is natural to introduce a coordinate that has these 
ramification points at and oo. So we define 

(4.1) x = t±l. 
v ' t-1 

Proposition 4.1. The Laplace transform of Dq 2(mi> M2) i s given by 

(4.2) F^{t x M)= E D^ 2 ( P i l ^ 2 )e-^ + ^ = -\og{l-z(x 1 )z{x 2 )) 

= log(ti - 1) + log(t 2 - 1) - log(-2(ti + t 2 )), 



where z(x) is the generating function of the Catalan numbers (3.1), and the variables 
t,w,x,z are related by (3.9), (3.13), and (4-D- We then have 

rD r. . \ . j rpD r. , \ dt\ • dt% dx\ ■ dx 2 dt\ ■ dt 2 



(4.3) W- 2 (t 1 ,t 2 ) = d 1 d 2 F^ 2 (t 1 ,t 2 )- {h _ t2)2 {xi _ x2)2 - {h + t2)2 . 
Proof. In terms of x = e w , the Laplace transform ( |4.2[ ) is given by 

(4.4) £ D 0l 2(/*i,A*2) e-^^+w^) 

1 ^ 1 /2j\ /2k\ J_ J_ ^ 1 /2 A /2fc\ _L 1_ 

" 4 v jJ U ; *y + jb jTin v jJ U ; xF 1 xip 1 

Since 

(4.5) dx = ^1 - cfe, 
we have 

(4.6) x- = ^±1 - = Z( " 2 ± 1} A. 

cix 1 — It z 2 — 1 cfe 

To make the computation simpler, let us introduce 

m=0 v 7 

This will also be used in Section [3 In terms of z and t we have 



(4.8) £ (z) 



1/ _^r^p m W 
2 \ dx J ^ m + l\m J x 2m ~ 



_ fii r 2m+l 

1 / z(z 2 + 1) d 



2\~ z 2 -l dz J z 2 -l At 

Note that 

/ d iUlA 1 /2A /2lfe\ 1 1 

- xi- V x 2 - 



dxi dx 2 J y4 ^ j + k\j J \ k J x 2 i x 2 2 k 
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+ 



1 



oo 
j,k=0 



X 



1 

2j+l 
1 



1 

„2k+l 



1 



(xi£ (xi) - l)(x 2 £o(x 2 ) - 1) + 2£o0ri)£o(>2) 



2z\Z2- 



1 + ZlZ 2 



+ 



z^zf + l) d 
z\ — \ dz\ 

In other words, we have a partial differential equation 

d 



■1) 

«a(*2 + 1) 



1 tfe 2 



(-log(l - ZiZ 2 )) . 



x\ 



for a holomorphic function in xi and £2 defined for |xi| ^2 and \x%\ >>> 2 . Since the first 
few terms of the Laurent expansions of — log (l — z(xx)z(x 2 )) using (3.7) agree with the 
first few terms of the sums of (4.4), we have the initial condition for the above differential 
equation. By the uniqueness of the solution to the Euler differential equation with the 
initial condition, we obtain (4.2). Equation (4.3) follows from differentiation of (4.2). □ 



In terms of the t-coordinate of (4.1), the Galois conjugate of t G S under the projection 
x : £ — >C is —t. Therefore, the recursion kernel for counting of dessins is given by 



(4.9) K D {tM) 



1 JTWi&Mi) _i 

2 W°(-t)-WR(t) 2 



1 



+ 



1 



t + ti t-ti 
1 



1 

64 



t + ti + t 



1 


1 


• dt 


t+1 t-1 




t-1 t+1 




1 \(* 2 - 


1)3 


1 






eft 



One of the first two stable cases (2.4) gives us 
= — I K D (t,t 



(4.10) Wii(t 



W^ 2 (t,-t) + 



dx • dx 



x 



Xl 



1 

2~7ri 



K D (t,h) 



dt ■ dt 



1 

128 



where the integration contour 7 consists of two concentric circles of a small radius and a 
large radius centered around t = 0, with the inner circle positively and the outer circle 
negatively oriented (Figure 4.1). The (g,n) = (0,3) case is given by 

W D 2 (t,t 1 )W D 2 (t,t2)W D 2(t,t 3 ) 



(4.11) W D 3 (t 1 ,t 2 ,t 3 ) 



1 

16 



1 

2~TTi 



dx(t) ■ dy(t) 
(t 2 -\) 2 {t-l)' 



dt 



(t+t 1 nt + t 2 ) 2 (t + t 3 y t 



dt\dt2dt 3 



1 

16 



f2 f2 4-2 

l l L 2 <-3 



dtidt 2 dtz 



Remark 4.2. The general formula (2.3) for (g, n) = (0, 3) also gives the same answer. This 
is because Wq 2 ac ts as the Cauchy differentiation kernel. 
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W^ 3 (h,t 2 ,t 3 ) = ^, J K D (t,h) W D 2 (t,t 2 )W D 2 (-t,t 3 ) + W D 2 (t,t 3 )W D 2 (-t,t 2 ) 

1 \(t 2 - l) 3 



1 

64 



-f —■ 



1 (t? - 1) 



32 tf 



t 2 V(* + *2) 2 (*-* 3 ) 2 (t + t 3 )*(t-t 2 y 

• dt\dt 2 dt 3 

1 1 

+ 



dt 



± d f t 3 (t 2 - l) 3 1 
16 5*3 V*l-*i 4 (t2 + t 3 ) 2 



I ( t 2 (t|-l) 3 1 



1 

j-2 j-2 j-2 

(- 1 t 2 tg 



i6dt 2 \tl-t\ t 2 {t 2 + hf 



dt\dt 2 dt 3 



1 

16 



dt\dt 2 dt 3 . 





/' 


-f— ^ 


TV' 






\ i 


.-y 




?-plane 



Figure 4.1. The integration contour 7. This contour encloses an annulus bounded 
by two concentric circles centered at the origin. The outer one has a large radius r > 
maxjgTv \tj\ and the negative orientation, and the inner one has an infinitesimally 
small radius with the positive orientation. 

Theorem 4.3. Let us define the Laplace transform of the number of Grothendieck's dessins 
by 

(4.12) F g D n (h, . . . = E D gtn (jj)e-^ + - + ^\ 



where the coordinate t{ is related to the Laplace conjugate coordinate Wj by 

tj - 1 tj + 1 



Then the differential forms 
(4.13) 

satisfy the Eynard-Orantin topological recursion 



W^ n (ti, ■ ■ ■ itfi) = df ■ ■ d n F® n (ti, . . . , t n ) 



(4.14) W° n (h,...,t n ) 



1 1 



1 1 

+ 



(t 2 - l) 3 1 



64 2iri J 1 \t + ti t - h J t 2 dt 
£ (Wg 2 (t, tj)W g , n -i{-t, t 2 ,...,tj,---,tn) + W£ 2 (-t, tj)W g>n -x(t, t 2 ,...,tj,---, tSj 
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stable 

+ w g \ n+1 (t,-t,t 2 ,...,t n ) + Yl w w,iji+i(*'^) w S,iji+iH.*j) 

91+92=9 
IUj={2,3,...,n} 

The last sum is restricted to the stable geometries. In other words, the partition should 
satisfies 2g\ — 1 + |/| > and 2g 2 — 1 + \ J\. The spectral curve £ of the Eynard-Orantin 
recursion is given by 

\x = z + \ 

\y = -z 

with the preferred coordinate t given by 

t -£±i. 

z - 1 

We give the proof of this theorem in the appendix. 

The problem of counting dessins is closely related to the counting problem of the lattice 
points of the moduli space M. g , n of smooth n-pointed algebraic curves of genus g studied 
in |59| 160] . Let us briefly recall the combinatorial model for the moduli space M g , n due 
to Thurston (see for example, [73]), Harer [36], Mumford [58], and Strebel [75], following 
|54l 155] . For a given ribbon graph T with e = e(T) edges, the space of metric ribbon 

graphs is M+ r ^/Aut(r), where the automorphism group acts by permutations of edges (see 
|544 Section 1]). When we consider ribbon graph automorphisms, we restrict ourselves that 
Aut(r) fixes each 2-cell of the cell-decomposition. We also require that every vertex of a 
ribbon graph has degree 3 or more. Using the canonical holomorphic coordinate system 
on a topological surface of \5i\ Section 4] corresponding to a metric ribbon graph, and the 
Strebel differentials [75] , we have an isomorphism of topological orbifolds [361 EE] 

(4.15) M g , n x M™ = R g>n 

for (g, n) in the stable range. Here 

TT < F) 

R w ~ 11 Aut(r) 

r ribbon graph 
of type (g,n) 

is an orbifold consisting of metric ribbon graphs of a given topological type (g,n). The 
gluing of orbi-cells is done by making the length of a non-loop edge tend to 0. The space 
R g ,n is a smooth orbifold (see [53 Section 3] and [73]). We denote by tt : R. 



g.n 



the natural projection via (4.15), which is the assignment of the perimeter length of each 
boundary to a given metric ribbon graph. 

Take a ribbon graph T. Since Aut(r) fixes every boundary component of T, they are 
labeled by [n] = {1,2 ... ,n}. For the moment let us give a label to each edge of T by an 
index set [e] = {1, 2, . . . , e}. The edge- face incidence matrix is defined by 

( 416 ) Av = Mie[n] tV e[e]> 

= the number of times edge n appears in face i. 

Thus 

a irj — 0; I, or 2, and the sum of the entries in each column is always 2. The T con- 
tribution of the space 7r~ 1 (/ii, . . . , fj, n ) = Rg )n (/j,) of metric ribbon graphs with a prescribed 
perimeter // = (/ii, . . . , fi n ) £ is the orbifold polytope 

{x G R% | ^ r x = fi} 

Aut(r) ' 
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where x = (li, . . . ,£ e ) is the collection of edge lengths of the metric ribbon graph T. We 
have 

(4.17) E ^ = E E °*tA = 2 £ v 

i£[n] ie[n] r;e[e] ?7G[e] 

Now let \x G be a vector consisting of positive integers. The lattice point counting 
function we consider is defined by 

|{xGZ" I A r x = /j}\ 



(4-18) N g>n (ji) = Yl 



|Aut(r)| 

r ribbon graph 
of type (g,n) 

for (c/,n) in the stable range ( fTO} l55j M\ IHO] ) . 

To find the spectral curve for lattice point counting, we need to identify the unstable 
moduli -Mo,i an d the ribbon graph space i?o,i- We recall that the orbifold isomorphism 



(4.15) holds for (g,n) in the stable range by defining R g , n as the space of metric ribbon 
graphs of type (g,n) without vertices of degrees 1 and 2. For (g,n) = (0,1), there are no 
ribbon graphs satisfying these conditions. Let Vj denote the number of degree j vertices in 
a ribbon graph T of type (g, n). Then we have 

j>i i>i 
where v is the total number of vertices of T. Hence 

(4.19) 2(2g -2 + n) = 2e-2v = - 2) Vj = - Vl + ^{j - 2) Vj . 

It follows that the number of degree 1 vertices v\ is positive when (g, n) = (0, 1). Thus we 
conclude that there is no spectral curve for this counting problem. 



Still we can consider the Laplace transform of the number (4.18) of lattice points of the 
moduli space M g ^ n with a prescribed perimeter length. We define for every stable (g, n) 

n 

(4-20) F g L n (t 1 ,... 1 t n )= £ NgMll-pr, 

Atez™ i=l Zi 

and the Eynard-Orantin differential forms by 

(4.21) W£ n (t u ...,tn) = d 1 ---d n F g L >n (*!,..., t n ). 

The following result is proved in |10| . with inspiration from [60J. 



Theorem 4.4 ( [10 ) . The differential forms W^ n {t\, . . . ,t n ) satisfy the Eynard-Orantin 



topological recursion with respect to the same spectral curve (3.13) and the recursion kernel 
\4-ty> starting with exactly the same first two stable cases 



(4-22) W L l{h) = _±_^_^f dtl; 



and 



1 



(4.23) W^ 3 (t u t2,t 3 ) = -_(! - -2-2-2 ) dhdt 2 dt 3 . 



16 V t{ t\ t 3 



Remark 4.5. It is somewhat surprising, because the spectral curve ( |3.13 ) has nothing to 



do with the lattice point counting problem. As we have mentioned, the [g, n) = (0, 1) and 
(0, 2) considerations for this problem do not produce the spectral curve. 
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Corollary 4.6. For every (g, n) with 2g — 2 + n > 0, we have the identity 

(4.24) W g D n (h, ...,t n ) = W g L n (h, ...,t n ). 

The differential form W^ n (t\, . . . , t n ) is a Laurent polynomial in t 2 , . . . , t n of degree 2(3g 
3 + n), with a reciprocity property 

(4.25) W£ n (l/h, l/t n ) = (-l) n tj ..-tl W g D n (t u ...,t n ). 

The numbers of dessins can be expressed in terms of the number of lattice points: 



(4.26) £> 9 ,„(/zi,...,/i n )= E"' Ell 



'> 2 1 1 



2^i \Xi I //< 

ii 



Hi 



N g , n (2h -m,--- ,2£ n -Hn)- 



Remark 4.7. The relation (4.26) appears in |62} Section 2.1] for an abstract setting. 

Proof. The Eynard-Orantin topological recursion uniquely determines the differential forms 
for all (g,n). Since W^(t) = W^ x (t) and W^ 3 (t 1 ,t 2l t 3 ) = W^ 3 (t 1 ,t 2 ,t 3 ), we conclude that 
W g D n (h, . . . , t n ) = W g L i(h, tn) for 2g - 2 '+ n > 0. 

By induction on 2g — 2 + n we can show that W? n (ti, . . . ,t n ) is a Laurent polynomial 
in t\, . . . ,t n . The statement is true for the initial cases (4.10) and (4.11). The integral 
transformation formula (4.14) is a residue calculation at t = and t = oo. By the induction 
hypothesis, the right-hand side of (4.14) becomes 

1 \ (t 2 - l) 3 1 



1 1 

64 2vri 

n 



1 



+ 



t + h t-ti 



t 2 



■dt' dh 



i=2 



(W£ 2 (t, t^Wg^-t, t 2 ,...,tj,...,t n ) + W£ 2 {~t, tj)Wg, n -l(t, t 2 ,...,tj,..., tn)^j 

stable 

+ w g D _ hn+1 (t,-t,t 2 ,...,t n )+ w£ iW+1 (t,t!)w£ M+1 (-t,tj) 

1 1 f (t 2 - l) 3 1 1 



91+92=9 
IUj={2,3,...,n} 



32 2ni J t 2 - t 2 t dt 



■ dU 



™2(t 2 + t 2 ) 

/ , 2 _ + 2\2 ^fl,n-l(*> *2j • • • j tj, . . ■ , t n ) dt ■ dtj 
3=2 [t V 



stable 



+ w g D _ lin+1 (t,t,t 2 ,...,t n )+ 



w. 



D 



91 



91+92=9 
JUj={2,3,...,n} 



Clearly the residues at t = and t = oo are Laurent polynomials in if , . . . , t^. 

Because of (4.24), we have 
(4.27) 

E °*»(aO]W^) = £ %Mf[rf(^) = (-ir E ^Hn< ! - 

M ez™ i=i ^ 1 ' ^ez™ i=i ^ * ' uz.z™ i=i 

where Xj = Zj + 1/zi- The Galois conjugation t— > — t corresponds to z^-l/z. Since 

w»(h, ...,t n ) = {-lyw^J-tu -t n ), 



the second equality of (4.27) follows. Take the residue of the left-hand side of (4.27) at 
Xi = oo for i = 1, . . . , n. On the right-hand side we take the residue at z% = for every i. 
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Then for every (/xi, . . . , fj, n ) G we have 
(4.28) 

1 

2tH 



Since 



|zi|=e J|z n |=e 



/'i 



n 

< n E 



1=1 



Zi + 



E 

^=0 



4: ' 1 



the residue of (4.28) comes from the term jUj — 2£i + vi = 0, and we have 



E E n(24-^)(^)^,n(24-//l,...,24-^n)- 



The reciprocity relation, and the degree of the Laurent polynomial, is the consequence 
of the following, which was established in [55J . 



Theorem 4.8 ([55]). The functions Fg n (ti, ... ,t n ) of (4-20) for the stable range2g—2+n > 
are uniquely determined by the following differential recursion formula from the initial 
values Fq 3 (tx , ti , £3) and (t\). 



(4.29) F ff yt l5 ...,i n ) 



1 

16 



ti 
-1 



E 



tj / (t 2 - 1) 3 d 



t 2 - 1 2 



F„ 



t 2 dt^ 9 ' 71 - 1 



{t,t[n)\{l,j}) 



t 2 at. 



F g,n-l( t [n]\{l}) 



+ 



+ E t 2 ( ,t' //■»-! 

1 (t 2 - 1 " 



(M[n]\{l,j}) 



J=2 

,2 1)3 ^2 



F g-l,n+l( u l, "2,i[ n ]\{l}) 



2 i 2 du\du2 

stable N 

31+92=9 / 
IUj=[n]\{l} 



dt. 



Here [n] = {1, 2, . . . , n} is an index set, and the last sum is taken over all partitions g\ +52 = 
g and set partitions I\J,J = [n] \ {1} subject to the stability conditions 2g\ — 1 + |/| > and 
2g2 — 1 + I J\ > 0. The initial values are given by 

1 (t + 1) 4 



(4.30) 
and 
(4.31) 



384 t 2 



t-A 



1 



F^ 3 (t 1 ,t 2 ,t 3 ) = -^(t l + l)(t 2 + l)(t 3 + l) fl + _L_ 



In the stable range Fg n (ti, . . . ,t n ) is a Laurent polynomial of degree 3(2n — 2 + n) and 
satisfies the reciprocity relation 

(4.32) F£ n (l/t 1 ,...,l/t n ) = F% n (t 1 , 
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The leading terms of Fg n {t\, . . . ,t n ) form a homogeneous polynomial of degree 3(2g — 2 + re), 
and is given by 

(4.33) ^( tl ,...,t n )^-g_L- £ (^•••r dn ) 9 , n n( 2 ^- 1 ) !! 

diH \-d„ j=l 

=3g—3+n 

where 

(r dl ---r dn ) g , n = ! i>i---^t 

JMg,n 

is the i^-class intersection number (see Section^for more detail about intersection numbers). 
The special value at ti = 1 gives 

(4.34) F g L n (l,l,...,l) = (-l) n X(M g , n ). 

This completes the proof of Corollary |4.6| □ 




5. The i/'-class interaction numbers on M g>n 



The crucial discovery of Konstevich [46] is the equality between the intersection numbers 
on the compact moduli space M.g t n and the Euclidean volume of the moduli space M.g jn 



of smooth curves using the isomorphism (4.15). The Feynman diagram expansion of the 



Kontsevich matrix integral relates the Euclidean volume with a r-function of the KdV 
equations. The Eyanrd-Orantin recursion for the V'-class intersection numbers is precisely 
the Dijkgraaf-Verlinde-Verlinde formula [15] of the intersection numbers. In this section we 
identify the spectral curve and the recursion kernel for the V'-class intersection numbers. 



As we have noted, the derivative of the recursion formula (4.29) is not the Eynard-Orantin 



recursion because the spectral curve is not defined by the unstable geometries. Indeed, we 
have i = 0. However, when we associate the number of lattice points with the V'-class 
intersection numbers on Mg t ni the unstable geometries do make sense. 
Let us recall a computation in |55} Section 4]. 



(5.1) N aMe 



(fj,,w) 



^ ^ |Aut(r)|^ x 

T ribbon graph /xSZ^ 1 V ' 1 
of type {g,n) 



6 



7 e(F) 



E 



r ribbon graph 
of type (g,n) 



|Aut(m ^ 



e (r) 



E 



r.mi 11 



|Aut(r)| 

r ribbon graph V edgef r) =l 

of type (g,n) of T 



E 



1 



Apx = /u}|e 



-{fj,,w) 



lAllKD II l_ e -(a v ,w) 
r ribbon graph v edge 

of type (g,n) of Y 



-(a v ,w) 



where is the incidence matrix of (4.16), a v is the 77-th column of Ai>, and (fi>,w) = 
HiWi + • • • + fJ-nWn- By comparing (4.20) and (5.1 ), we see that we are substituting e Wi = Zi 
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in this computation. Therefore, we obtain 



(5.2) 



1 



n 



|Aut(r)| J- J- rp» z air > - 1' 

r ribbon graph r\ edge 1 Ll — 1 1 

of type (g,n) of Y 



Thus the series (4.20) in Zi converges for \zi\ > 1. Since Z{ = fj^zr, the tj— )-oo limit picks 
up the limit of (4.20) as z~ >1, and hence the information of N gjn (/j,) as /x— >oo. Since the 
orbifold isomorphism (4.15) is scale invariant under the action of M+, making the perimeter 



length /i large is the same as making the mesh small in the lattice point counting. Hence at 
the limit we obtain the Euclidean volume of M g , n considered by Kontsevich in |46j . This is 
why we expect that (4.33) holds. Let us now consider the limit of the spectral curve (3.13) 



as t— >oo. First we have 



z + 



1 



V 



z 

-z = -1 



2 + 



t 2 - 1 
2 



t- 1 



Ignoring the constant shifts of x and y, we obtain for a large t 
(5.3) 



x 

y 



2 
t ■ 



Hence the spectral curve is given by the equation x = y 2 . We use t as the preferred 
coordinate. 

We now compare the Eynard-Orantin recursion with respect to this spectral curve and 
the Witten-Kontsevich theory. We use ( |4.33 | and define 

^g,n(*l) ■ • ■ i*n) = d\ ■ ■ ■ d n F gn (ti, . . . , t n ) 

2dj 



■1 



(5.4) 



E (Td 1 ---T dn ) gtn l[(2d j + l)\ 



22g-2+n * ' ™i ' <*ni g,n X J. \ J ' ~'" \ 2 

diH \-d n j=l 

=3g— 3+n 



(-l) r 



1 • 



lg2g-2+n 3,™ v 

where w g ^ n (t\, . . . ,t n ) is the coefficient of the Eynard-Orantin differential form normalized 

by the constant factor ^g^+n ■ Note that wjf n (tx, . . . , t n ) is a polynomial in t 2, s with 
positive rational coefficients for (g,n) in the stable range. For (g,n) = (0, 1) and (0,2), we 
have 



(5.5) 
(5.6) 

Therefore, 
(5.7) 



( T k)o,_ 



Jk+2,0 



(T kl T k2 )o,2 = (-l) kl : 



h + k 2 = -1. 



■1 



W&(t) = I ^ T (r_ 2 )(-3)!! t-'dt = ^dt = ydx, 



16 



in agreement with the spectral curve x = y (5.3). Similarly, we have 

2d+l , ± \ -2d-l 



(5.8) if 2 (i!,i 2 ) = f>l) d (2d " 1)H(-M - 3)!! f |) 
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t—(*) 



2d+l 



log 1 



log 1 



and hence 
(5.9) 



wF 2 (h,t 2 ) 



dt\ ■ dt 2 1 dx\ ■ dx2 



(ti - i 2 ) 2 2 (xi 
As a consequence, the recursion kernel is given by 



X2 



\2- 



(5.10) 



K K (t,h) 



1 



2 W + ii 



+ 



i - U 



32 dt 



dU 



since /^''^^ does not contribute to the kernel. The Eynard-Orantin recursion for the 



Euclidean volume then becomes 

K \n, ■ ■ ■ , <>n) 

1 



(5.11) W^ n (h,...,t r 



1 

2tH 



loo 



+ 



1 



t + ti t-hj 64 dt 



t 4 1 



dti 



w« 1>n+1 (t, 



-t, t 2 , ■ ■ ■ , t r 



dt ■ dti 




-l( t, t 2 , ■ ■ ■ , tj, . . . , t n ) 



3,n-l 



(M2, 



tj, ■ ■ ■ , t n ) 



stable 



+ E w ?,[/i+i(*.*0wJ,| JI+1 H,tj) 

9i+ff2=g 

/UJ={2,...,7J.} 

where the integral is taken with respect to a large negati vely oriented circle 7oo that encloses 



any of ±ti, . . . , =ti n . This is the larger circle of Figure 



4.1 



Here again dx ^ dx z does not 

° (XI— X2) 



contribute in the formula. Since the coefficients Wg n (t\, . . . ,t n ) in the stable range are 
polynomials, the poles of the integrand of (5.11 ) in the integration coutour are at t = ±ij's. 
Therefore, we can perform the integral in terms of the residue calculus at poles t = ±ij. 
First let us get rid of the factor 1/16 29 ~ 2+ ™ from (5.11). Since the recursion is an induction 
on 2g— 2+n, we have an overall factor 16 adjustment on the right-hand side. The integration 
contour is negatively oriented, so the residue calculation at t = ±tj receives uni versa lly the 
negative sign. This sign is exactly cancelled by the choice of the sign of w^ n in (5.4). Thus 
the result of residue evaluation of (5.11) is 



(5.12) < n (ti,. 



i tn) 



hw g _ ln+1 (ti,ti,t 2 , 



■ ) t n ) 



stable 



l \ E <,iii+i(*i»*')<,ui+i(<1'*j) 



91+92=9 
IUj={2,...,n} 



tj + t 2 



J=2 



Hi 



3 ...K 
fiy w g,n-l 



(t\ , . . . , tj , . . . , t n ) 
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1 " 



3=2 







9 




- 


H 


" at 


*=-tJ ( 



^ W g,n—l(ti t%t ■ ■ ■ ) • • • > ^n) 



stable 



w£l,„+l(tl,tl,*2, •••,*»») + J] ^,|/|+l(*l^/)^£,|J| + l(*l 5 *j) 



+ E 

3=2 



0fc 



91+92=3 
7UJ={2,...,n} 



+2 _ f 2 



4„..K 



This is the same as [TOj, Theorem 5.2]. 
Let us adopt the normalized notation 

n 

(5-13) (a dl ■ ■ ■ a dn ) g>n = (r dl ■ ■ ■ r dn )^ n J](2^ + 1)!! 

i=l 

to make the formula shorter. Then 

(5-14) w£ n (h,...,t n ) = (< T d 1 ---<Td n ) 9tn 

di,...,d„ j=l 

The DVV formula |15j for the Virasoro constraint condition on the ^-class intersection 
numbers on M. g , n reads 

n ^ n 

(5.15) (<T fe Yl<Tdi)g, n = x Yl ( <r «°*II '*^-i,™+i 



i=2 



a+fe=fc-2 
stable 



i=2 



92,|J|+1 



a+fe=fc-2 91+92=9 

7UJ={2,...,n} 



+ ^(2d J + l)(^ +dj ._ 1 [] ^), iJt _ 1 . 

3=2 

We thus recover the discovery of |25j : 

Theorem 5.1. The Eynard-Orantin recursion formula for the spectral curve x = y 2 is the 
Dijkgraaf-Verlinde-Verlinde formula [15] for the intersection numbers (r dl ■ ■ ■ T dn ) g n on the 
moduli space Mg, n of pointed stable curves. 
Proof. We extract the coefficient of 



(5.16) 



tf n t? j 

3=2 



in (5.12) and compare the result with (5.15). It is obvious that the fifth line of (5.12) 
produces the first and second lines of (5.15). 



To compare the last lines of (5.12) and (5.15), we consider the case \tj\ < \ti\ for all j > 2 
in (5.12). We then have the expansion 
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The (5.16)-term of the last line of ( |5,12 ) has two contributions. The first one comes from 



d_ 



ffij ^ ( l| ) w ^n-l( t l,t2,---,t j ,...,t n ) 



m=0 



Since w^ n _ l {t\,t2, ■ . . , tj, . . . , t n ) does not contain tj, we set m = dj to produce the right 
power 2dj of tj. The power of t\ has to be 2k. Thus from w^ n _ 1 we take the term of 
+2 j 2 , whose coefficient is ((Tk+dj-i Yli^i j a d z )- The total contribution from the first 



kind comes from the differentiation, which gives 2m + 1 = 2dj + 1 




The second possible contribution for the (5.16)-term may come from 

±5 00 



However, this term does not produce t\ , and hence does not contribute to the (5.16)-term. 



This completes the proof of Theorem 5.1 □ 



6. Single Hurwitz numbers 

What is the mirror dual of the number of trees? The answer we wish to present in this 
section is that it is the Lambert curve. This analytic curve serves as the spectral curve for the 
Hurwitz counting problem, and comes up from the the unstable geometries {g, n) = (0, 1) 
and (0, 2) via Laplace transform. 

A Hurwitz cover is a holomorphic mapping / : C — > P 1 from a connected nonsingu- 
lar projective algebraic curve C of genus g to the projective line P 1 with only simple 
ramifications except for 00 £ P 1 . Such a cover is further refined by specifying its pro- 
file, which is a partition \x = (jx\ > fj,% > • • • > /i„ > 0) of the degree of the covering 
d = \fJ,\ = fj,i H — ■ + [i. n - The length £(fi) = n of this partition is the number of points in the 
inverse image / _1 (oo) = {p\, . . . ,p n } of 00. Each part /ij gives a local description of the 
map /, which is given by u 1 — > in terms of a local coordinate u of C around pi. The 
number hg jfl of the topological types of Hurwitz covers of a given genus g and a profile fi, 
counted with the weight factor l/|Aut/|, is the single Hurwitz number we shall deal with 
in this section. 

The deformations of a Hurwitz cover / are obtained by moving the branch points (i.e., 
the image of the ramification points) on P 1 \ {00}. Thus h g ^ counts the number of Hurwitz 
covers with prescribed (i.e., fixed) and labeled branch points. On the other hand, the 
preimages of 00 on C are labeled only by the parts of \x. Therefore, a more natural count 
of Hurwitz cover is 

la -n tj f \ |Aut(/x)| 

(fU) = (2,-2 + n + H)! • k ^ 

Here, 

(6.2) r = r{g,^) A =2g-2 + n+\n\ 

is the number of simple ramification points of / by the Riemann-Hurwitz formula, and 
Aut(/i) is the group of permutations of equal parts of the partition /j,. 

One reason that explains why single Hurwitz numbers are interesting is a remarkable 
formula due to Ekedahl, Lando, Shapiro and Vainshtein [20, 33, 48, 65J that relates Hurwitz 
numbers and Gromov-Witten invariants. For genus g > and a partition [i of length 
i{fj) = n subject to the stability condition 2g — 2 + n > 0, the ELSV formula states that 
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Mr A v 



(6.3) *M = U%L rr 



A^(l) 



[1 - Liiifji) 



9 n m+ki 

i=o fcx,...,fc n >o i=i 

where M g , n is the Deligne-Mumford moduli stack of stable algebraic curves of genus g with 
n distinct smooth marked points, (1) = 1 — ci(E) + - • • + (— l) 9 c 9 (E) is the alternating sum 
of the Chern classes of the Hodge bundle E on M g , n , ipi is the i-th tautological cotangent 
class, and 

(6-4) (r kl • • • r kn c,(E)) = [ ^ ■ ■ ■ ipg n Cj(E) 

J Mg,n 

is the linear Hodge integral, which is unless k± H + k n + j = 3g — 3 + n. 

The Deligne-Mumford stack M g , n is defined as the moduli space of stable curves satisfying 
the stability condition 2 — 2g — n < 0. However, single Hurwitz numbers are well defined 
for unstable geometries (g,n) = (0, 1) and (0,2), and their values are 

jd— 3 d^^ 1 
(6.5) Ho ((d)) = — — = — and H ((^,^)) 



{d - 1)1 d\ ' ' m + H2 H\\ /i 2 ! 

The ELSV formula remains valid for unstable cases by defining 

A V (1) _ 1 



(6.6) 
(6.7) 



Mo i 1 - # d 2 ' 



A V (1) 1 



Mo,2 



(1 - /Ul^l)(l - ^2^2) Ml+/^2 



Let us examine the (<?, n) = (0, 1) case. We wish to count the number of Hurwitz covers 
/ : P 1 — 7-P 1 of degree d with profile (i = (d). If d = 2, then /(u) = u 2 is the only map, 
since r = 1 and the two ramification points can be placed at u = and n = 00. The 
automorphism of this map is Z/2Z. We now consider the case when d > 3. First we label 
all branch points. One is 00, so let us place all others, the images of simple ramification 
points, at the r-th roots of unity. Here r = d — 1. We label these points with indices 
[r] = {1,2, ... ,r}. Connect each r-th root of unity with the origin by a straight line (see 



Figure 6.1). Let * denote this star-like shape, which has one vertex at the center and r 
half-edges. Then the inverse image / _1 (*) is a tree-like shape with d vertices and rd half- 
edges. Here we call each inverse image of a vertex of / _1 (*)- If / is simply ramified at p, 
then two half-edges are connected at p and form a real edge that is incident to two vertices. 
Since f{p) is one of the r-th root of unity, we give the same label to p. Thus all simple 
ramification points are labeled with the index set [r]. Now we remove all half-edges from 
/ _1 (*) that are not made into an edge, and denote it by T. It is a tree on P 1 that has d 
vertices and r = d—1 edges. Note that except for the case d = 2, the edge labeling gives a 
labeling of vertices. For example, if a vertex x is incident to edges i\ < 12 < • • • < ik, then 
x is labeled by i±i2 ■■■ik- 

Conversely, suppose we are given a tree with d labeled vertices by the index set [d] = 
{1, 2, . . . , d} and r = d—1 edges. At each vertex we can give a cyclic order to incident edges 
by aligning them in the increasing order of the labels of the other ends of the edges. Thus 
the tree becomes a ribbon graph (see Section [3]), and hence it can be placed on P 1 . Then 
by choosing the midpoint of each edge as a simple ramification point and each vertex as a 
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Figure 6.1. Counting the genus single Hurwitz numbers with the total ramifi- 
cation at oo. 



zero of /, we can construct a Hurwitz cover. Recall that the number of trees with d labeled 
vertices is d d ~ 2 . Therefore, 

d d - 2 

is the number of trees with d unlabeled vertices. 

Fix an n > 1 , and consider a partition fi of length n as an n-dimensional vector 

H = (//i, ...,n n ) e z+ 

consisting of positive integers. The Laplace transform of H g (fj,) as a function in /j,, 

(6.8) Hg tn (wi, . . . , w n ) = H g (n)e 

is the function we wish to compute. Note that the automorphism group Au t(//) acts trivially 
on the function e~ ^^n(w n +i)) ^ w j 1 j c ] 1 explains its appearance in (6.1 ). The reason 
for shifting the variables Wi i — > Wi + 1 is due to the asymptotic behavior 

At! V2tt 

as /i approaches to oo. This asymptotics also suggests that the holomorphic function 
H g>n (wi, . . . , w n ) is actually defined on a double-sheeted overing on the uij-plane, since 
y/vJi behaves better as a holomorphic coordinate. 

Following [231 [57], we introduce a series of polynomials £ n (t) of degree 2n + 1 in t for 
n > by the recursion formula 

(6-9) £n(t)=* 2 (*-l)^n-l(*) 

with the initial condition £o(t) = t — 1. This differential operator appears in |31j . The 
functions £-i(t) and ^o(^) also appear as the two fundamental functions in [84J that generate 
his algebra A. These polynomials are introduced to make the computation of the Laplace 
transform (6.8) easier. 
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Proposition 6.1 ( |121 124j ). Let us introduce new coordinates 



(6.10) 



x = e 



OO ,,1 

^ 1 



11=1 



Then the inverse function of z = z(x) is given by 

(6.11) x = ze 1 ~ z , 
and the variables z and t are related by 

t - 1 

(6.12) z ■■ 
Moreover, we have 



t 



(6.13) 



for n > 0. 



^) = E^r e ^ +1) = E 



fi+n 



//! 



Proof. The infinite series (6.13) has the radius of convergence 1, and for |x| < 1, we can 
apply the Lagrange inversion formula to obtain ( 6.11| ). Since the application of 

d d 9 . d 

— r~ = x — = t 2 (t - 1 — 
dw ax dt 



n-times to X^^Li Tjf e ^ w+1 ^ produces Yl^=i e we obtain (6.9). If we extend 

(6.13) formally to n = —1, then we have z = £_i(i). To obtain the expression of z as a 
function of t, we need to solve the differential equation 



t\t-r 



d_ 
dt 



z = t-l. 



Its solution is z = c — \. Since x = <s=> z = and x = ==>• i = 1, we conclude that the 
constant of integration is c = 1. Thus z = 1 — 1/i. □ 

Remark 6.2. The relation between our z as a function in x and the classical Lambert 
W- function (see for example, [12]) is 

z(x) = —W(—x/e). 



Because of the ELSV formula (6.1 ), the Laplace transform of H g (n) becomes a polynomial 
in ti, . . . , t n for (g, n) in the stable range. The result is 



(6.14) F^ih, . . . ,t n ) = H g , n {w(t l ),...,w(t n )) 

= £ H g [n)e-^+ l 



E 



£ (T kl ...r kn A^(l)) Jja_^ e -0i(«i+i)+...+^(«.+i)) 

/UEZY&iH hfc„<3g-3+n 1=1 



e n4(i»). 

felH hfcn<3g— 3+n i=l 

The Laplace transform ( 6.14 ) is no longer a polynomial for the unstable geometries (g, n) ■ 
(0, 1) and (0, 2). Wei use (16.51) to calculate Fq X and F^ 2 . 
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Theorem 6.3. The Laplace transform of the unstable cases (g,n) = (0,1) and (0,2) are 
given by 



(6.15) 
and 
(6.16) 



F&ihM =log 



^1 - ^2 



- x 2/ 

where ts related by \6.11]j and §6.1 Sty . 

Proof. The (0, 1) case is a straightforward computation. 



(zi + z 2 ) + 1, 



F H i(t) = E e " v = E e " v = e "- 2(t) 



fc=cZ d=l 

This is a solution to the differential equation 

d 



t 2 (t-l)-^ 2 (t) = ^ 1 (t) = z 



t-1 



Therefore, £-2(i) = c — \ h for a constant of integration c. Here again we note 

t = l^z = 0^x = 0^ |_ 2 (t) = 0. 
This determines that c=\. Thus we have established (6.15). 



Since 



Ml A 4 2 



A t l:A 1 2>l 



fJ,\+fJ,2 £H ! 



and since z = £_i(t), (6.16) is equivalent to 
(6-17) E 



1 /; M2 



Ml,M2>0 

(Mi,M2)^(0,0) 



Hi + /i 2 Mi 



^2! 



■x 1 x 2 



log e E 



-k x l 



37 1 



fe=l 



fc! 



xi - x 2 



where \x\\ < 1, |x2| < 1, and < \x\ — x 2 \ < 1 so that the formula is an equation of 
holomorphic functions in x\ and x%. Define 



{xx,x 2 ) 



clef 



1 



,/'2 



Ml _-Ati . M2_ -jta . ^1^2 



/ii,/i 2 >0 
(Mi,/x 2 )^(0,0) 



/ii + ^2 Ml' 



M2! 



log E 



00 jjfc-l 



1-fc x l 



3y o 



\k=l 



k\ 



x% - x 2 



Then 



•°)= E ^Ve-^x«-log(E 

/Ltl>l \fc=l 



e fc • x fc 1 



|_i(t) - log ( j - 1 = 1 - - - log (l - - ] + loo ,. _ 1 



1 



log 1 



1 



w = 
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7 

Here t is restricted on the domain Re(t) > 1. Since 



x = e ~ w = ze 1 " 2 = I 1 — - ) e*. 



00 hk-l 



x\ — — log e y - ; e 



-fc ^1 



k=l 



t\{tl - 1)^- log (tl) " C-l(*2)) " XI— log(xi - X 2 ) 

2 lX 9 , / 1 , 1 \ X\ 



dti \ h t 2 J xi - x 2 

ht 2 (ti - 1) xi 



h — t 2 xi - x 2 



we have 



d d \ . Afc 



\ ^ -fe 



x i o 1" x 2t^ — log e } — — — e 



<9xi 9X2/ \ ^— ' fc! Xi — x 2 

\ k=l / 

M2(*i - 1) - ^1^2(^2 - 1) _ gl - %2 
il - *2 a;i — x 2 

= ht 2 -i = io(h)i (t 2 ) + e (*i) + io(t 2 ). 

On the other hand, we also have 

° , „_ \ 1 ^1 r -Ml ^2 -M2 r m r M2 



dxi dx 2 y ^— ' Ml + A*2 Ml' M2! 

A t l,M2>0 

(Ml, ^2)^(0,0) 



y M1 + M2 _ mT e - w . /ff e -M2 . x mi x M2 

(^1^2)^(0,0) 



&(tl)6(t2)+eo(*l)+&(*2). 



Therefore, 



(6.18) (xi^ + x 2] ^y { x u x 2 ) = 0. 

Note that c/>(xi,X2) is a holomorphic function in xi and x 2 . Therefore, it has a series 
expansion in homogeneous polynomials around (0,0). Since a homogeneous polynomial in 
xi and X2 of degree n is an eigenvector of the differential operator xigf- + ^2gf~ belonging 
to the eigenvalue n, the only holomorphic solution to the Euler differential equation (6.18) 
is a constant. But since <f>(xi,0) = 0, we conclude that <^(xi,X2) = 0. This completes the 
proof of (6.17), and hence the proposition. □ 

Definition 6.4. We define the symmetric differential forms for all g > and n > by 

(6.19) W* n {ti,...,t n ) = di---d n F^ n (ti,...,t n ), 
and call them the Hurwitz differential forms. 
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The unstable cases are given by 
(6.20) W^{h) = diF,fi(ii) = 73 dh = - dx, 



q x 



and 



(6.21) W&fah) = hfaFgJtutv) = did 2 [log (21 - z 2 ) - log(xi - x 2 )\ 



d\d 2 



did 2 [log (ti - t 2 ) - log(xi - x 2 )] 

dti ■ dt 2 dx\ ■ dx2 



~ (ti-t 2 f { Xl -x 2 f 

We note that all quantities are expressible in terms of z, or equivalently, in t. Now Defini- 
tion 2.1 tells us that the spectral curve £ of the single Hurwitz number is 

I x = ze x ~ z 

6 - 22 \ Z Z -l 

[y= I = e • 

The Lambert curve £ defined by x = ze 1 ^ z , which is obtained by the Laplace transform 
of the number of trees, is an analytic curve and its x-projection has a simple ramification 
point at z = 1, since 

dx = (1 - z)e l ~ z dz. 

The t-coordinate brings this ramification point to t = oo. Let z (resp. t) denote the unique 
local Galois conjugate of z (reps. t). We also use 

(6.23) t = s(t), 
which is defined by the functional equation 

(6.24) ( 1 - 1 ) e\ = ( 1 - ) eW5. 



tj V s (*), 

Although the Galois conjugate is only locally defined near the branched point t = oo, we 
consider s(t) as a global holomorphic function via analytic continuation. For Re(t) > 1, 
(6.24) implies 

w(t) 



io * = -f--V--Uv- 

° gX ~ \ t n t n } ^ t n ' 

\ n=l / n=2 



When considered as a functional equation, ( 6.24[ ) has exactly two solutions: t and 

< 6 - 25 > »w = - t n + il 4 - 2 + ^ 3 + ^ 4 + -- 

This is the deck-transformation of the projection n : £ — )■ C near t = oo and satisfies the 
involution equation s(s(i)) = t. It is analytic on C \ [0, 1] and has logarithmic singularities 
at and 1. 

Let us calculate the recursion kernel. Since 

dx 1 — z , dt s'(t)dt 

— = dz — — 



x : t 2 (t-l) s(i) 2 (s(i)-l)' 

we have 

(*M K»(+ + \ 1 I t S(t) W H 2 (;h) If 1 1 \ t\t-l) 1 

(6 - 26) * (t ' * l} = 2 - Wb,i(*) = 2 " T^fe ' dt ■ dh 
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1 



1 



1 



2 \t-ti s(t) -tij s(t) - t 



ts(t) t\t - 1) 



dt 



■ dh. 



Theorem 6.5 ( |244 I57j). The Hurwitz differential forms (6.19) for 2g — 2 + n > satisfy 
the Eynard-Orantin recursion: 



(6.27) W g H n (t 1 ,...,t n ) = ±.f K H (t,h) 



W g \ n+l (t,s(t),t 2 ,...,t n ) 



No (0, l)-terms 

91+92=9 
7UJ={2,...,n} 



where joo is a negatively oriented circle around oo whose radius is larger than any of \tj \ 's 
and \s(tj) \ 's. 



Remark 6.6. The recursion formula (6.27) was first conjectured by Bouchard and Marino 
in [7j. Its proofs appear in [H [241 157] . The method of [3] is to use a matrix integral 
expression of the single Hurwitz numbers. The idea of [21, 57J is that the Laplace transform 
of the cut-and-join equation of [ 291 EZj is the Eynard-Orantin recursion. The cut-and-join 
equation takes the following form: 

(6.28) r(g,v)H g (n) = + fJ,j)H g (fj,(i, j), m + fij) 

i<j 

( 



i=l a+/3=m 



\ 



Hg-i{^(i),a,P) + ^2 H gi{vi,oi)Hg 2 {v2,P) 



V 



9l+92=9_ 



I 



Here p is a partition of length n, and fi(i) and /J-(i,j) indicate the partition obtained by 
deleting parts of \i. 

Remark 6.7. As we have seen above, Hurwitz numbers in ( |6.5| ) determine the shape of the 
recursion formula (6.27). Since the recursion gives the Hurwitz numbers for all (g,n), we 
have thus established that unstable (g,n) = (0,1) and (0,2) Hurwitz numbers determine 
all other single Hurwitz numbers. 



It is important to check if the formulas (2.4) and (2.6) agree with the geometry. From 
the definition (6.14) we calculate 

F$ 3 (t h h,t 3 ) = (T T To} Jo(tl)Ut2)Ut3) = (h ~ 1)(*2 " 1)(*3 " 1), 

which yields 

(6.29) W£ 3 (h,t 2 ,t 3 ) = dhdt 2 dt 3 . 

Since 



dx{z) ■ dy(z) = (1 — z)dz ■ dz 



dt ■ dt 



from (6.22) and (6.11), the general formula (2.6) yields 



W? 3 (t u t 2 ,t 3 ) 



1 

"2vri 



Too 



W^jtMW^jtMW^tM) 
dx{t) ■ dy{t) 
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1 



r(lt 



dtxdt^dt^ = dtxdt-idtz 



(t - hnt - t 2 y(t - t 3 f 

in agreement with geometry. Here we calculate the residue at t = oo. Although 

t ji > ■ ■ n dt ■ dti dx ■ dxi 



W»Jt,ti 



{t-Uf (x-xi) 21 



the second term does not contribute to the integral. This is because as t— >oo, we have x— >1, 
and dx ■ dxij{x — xi) 2 has no pole at x = 1. 
Similarly, 



1 



F&ih) = <n)i,iCi(ti) - (roXd^oih) = - - 1), 

and thus we have 
(6.30) 



Wg(t 1 ) = ^(*i-l)(3t 1 + l)dti. 



On the other hand, the general formula (2.4) gives 

dx(u) ■ dx{v) 



WK(u,v) + 



x(u) — x(v)) 2 _ 

dt ■ s'(t)dt 



u=t 
v=s(t) 



l 

2vri 



1 / 1 _ 1 \ ta(t) 2 , _ 

2 V* — *1 s{t) - h J s(t) - t [ \t-s(t)) 2 

tis(h) 2 







2vri 


A[o,i] 2 ^ 



l/l -s(tl)) 3 

1 



s(ti) 2 (s(ti) - 1) dtl 



1 



ts(f)_ 2 s'{t)dt 

<t)-t [t ij (i-s(t)) 2 



where 7[o,i] is a contour circling around the slit [0, 1] in the t-plane in the positive direction. 




Figure 6.2. The contours of integration. The outer loop 700 is the circle of a large 
radius oriented clock wise, and 7[o,i] is the thin loop surrounding the closed interval 
[0, 1] in the positive direction. 
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Note that the integrand of the last integral is a holomorphic function in t on 7[o,i], hence 
it has a finite value. It is also clear that as t\— >oo, this integral tends to 0, because 7[o,i] is 
a compact space. Therefore, we conclude that 



w£(ti) 



tis(ti 



(ii-^i)) 3 



s(ti) 2 (s(ti)-l) dti + 0(l/h) 



1 

12 



1 

24 



dti + 0(l/h 



since s(t) = — t + 2/3 + 0(l/t 2 ). It agrees with (6.30) because of the following 
Lemma 6.8. A solution to the topological recursion (6.21) is a polynomial in t\. 



Proof. The independence of Wg n {t\, ■ . . , t n ) only comes from the factor 
1 1 



t - ii s{t) - h 



1 1 1 

-t ++ 3* 



2 2\ 1 
3 h + 9)¥ 



+ \4 



2 22 \ 1 
3 tl + 135 J ^ + 



in the recursion kernel (6.26). Since each coefficient of the ^-expansion of K H (t,t\) is a 
polynomial in t\, Lemma follows. □ 



7. The stationary Gromov-Witten invariants of P 1 

In this section we study the generating functions of stationary Gromov-Witten invariants 
of P . The conjectural relation between these invariants and the Eynard-Orantin topological 
recursion was first formulated in [62] . We identify the spectral curve and the recursion kernel 
using the unstable geometries. 

Morally speaking, the space P 1 we are considering here appears as the zero section of a 
Calabi-Yau threefold known as the resolved conifold X, which is the total space of the rank 2 
vector bundle P i (— 1) © Opi (— 1) over P . Let L C X be a special Lagrangian submanifold 
|39} 142] . Then the intersection L n P 1 of the special Lagrangian and the zero section is a 
circle on P . If we holomorphically embed a bordered Riemann surface with n boundary 
components into X in a way that each boundary is mapped to a distinct circle on P , then 
the whole Riemann surface is necessarily mapped to P . Thus we are considering open 
Gromov-Witten invariants of P . And if we make these circles on P 1 small and centered 
around n distinct points of P , then we are naturally led to the stationary Gromov-Witten 
invariants of P 1 . 

So our main object of this section is the Laplace transform of the stationary Gromov- 
Witten invariants 

OO 71 Th 

(7-1) F^ n ( Xl ,...,x n )= £ (r Ml H---r Mn H) 9in n^II^ 

jitl,...,ju n =0 i=l i=l 

where oj G ^-oClP 1 ) is the point class generator, and 

(7.2) <vH-^« n = I Vr^M-.-Vre^M 

J[M B>n ^,d)]^ 

is a stationary Gromov-Witten invariant of P . More precisely, JA g>n (¥ ,d) is the moduli 
stack of stable morphisms from a connected n-pointed curve (C,p\, . . . ,p n ) into P 1 of degree 
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d such that f(pi), i = 1, . . . , n, are distinct, and evi is the natural evaluation morphism 
e Vi : M g , n (F\d) 3 [/, (C, Pl , . . . ,p n )] .— > f( Pi ) E P 1 . 



The Gromov-Witten invariant (7.2) vanishes unless 

(7.3) 2g - 2 + 2d = m + • • • + ^. 



The sum in (7.1) is the Laplace transform if we identity 
(7.4) 



x = e 



The extra numerical factor n^A'j! * s included in (7.1) because of the polynomial growth 
order of 



(7.5) 



(r w (w)---r /ln (w)) fl)n JJ/z i 



i=l 



for large \x that is established in [66J. Indeed (7.5) is essentially a special type of Hurwitz 
numbers that counts the number of certain coverings of P . 

To determine the spectral curve and the annulus amplitude, we need to consider unstable 
geometries (g, n) = (0,1) and (0,2). From [66] we learn 

/i \ 2 

(7-6) (t>i( w )>o,i = <7Jw-2(w))o,i = 



To compute a closed formula for 



<(*) = E (^H)o,i w'^ttt = E^p 1 ^ 



/xi=Q 



we notice that the generating function of Catalan numbers (3.5) 

oo ^ 
z (^) = E Cm ~2m+l 



m=0 



provides again an effective tool. Thus we have 



(7.7) x 



n oo 

l)< 1 (x) = -2E 

' d=l 



(2d-2)\ 1 



{d-l)\d\ x™- 1 



- 2 E 

m=0 



(2m)! 



1 



m 



+ l)!m! x 2m+1 



-2z(x). 



The advantage of using the Catalan series z(x) is that we know its inverse function (3.8). 
Using (4.6), we see that (7.7) is equivalent to 

' z 3 + z d 



(7.8) 



z 2 — 1 dz 



l)F? 1 {z) = -2z. 



The solution of (7.8) is given by 



2 



z + - J log(l + z 2 ) + c f 2 + - 



with a constant of integration c. Since 

z-^0 x— >oo Fq\^0, 
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we conclude that c = 2. We thus obtain 
(7.9) 

and therefore, 
(7.10) 



F^(z) = 2z- (z+^j log(l + z 2 ), 



o,i i 



log(l + z 2 )d{ z + 



Theorem 7.1. The spectral curve for the stationary Gromov-Witten invariants of P 1 is 
given by 

l 



(7.11) 



x = z + ^ 
= -log(l + z 2 ). 

Remark 7.2. Since = has two zeros at z = ±1, we also use as our preferred coordinate 

/„ x z + l t + 1 

(7.12 t = <^z = . 

v ' z-1 t-1 

The lo g sin gularity on the i-plane is the right semicircle of radius 1 connecting i to —i (see 
Figure 7.1). The expression of Wq\ in terms of the preferred coordinate is 

'2(t 2 + l) 



(7.13) 



<(*) 



SI 



(t 2 



1) 



log 



(t-1) 



dt. 




Figure 7.1. The spectral curve for the stationary Gromov-Witten invariants of P 1 
is the complex i-plane minus the semicircle. 



Remark 7.3. The function x = z + - is expected here, since it is the Landau-Ginzburg 
model that is homologically mirror dual to P 1 [2]. 

Remark 7.4. The Galois conjugate of x = z + - is globally defined, and is given by 
(7.14) t\ — >t = -t. 
Remark 7.5. Since 



oo 



v ; i-z(x) ^ y ' ^ v m \ x n+i 

y ' k=0 n=0 VL 2J/ 
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we can express t in the branch near t = — 1 as a function in x. The result is 

z(x) + l _ „ 2 |W n \ 1 

n=0 ^f. 

which is also absolutely convergent for |x| > 2. 



(7-16) t + l = ^±i + l = 2- - -. =-V2(,„,, 

v ; z(x)-i i-^(x) VLijy^ n+i 



Remark 7.6. We are using the normalized Gromov-Witten invariants (7.5) to compute 
the Laplace transform (7.1). If we did not include the fj,\ factor in our computation of the 
spectral curve, then we would have encountered with the modified Bessel function 

oo 1 

W = E T^r 



oo i 

(m! 

m=l 



|\2~ ' 



instead of z(x), in computing (7.9). We note that Iq(2x) appears in [18] in the exact same 
context of computing the Gromov-Witten invariants of P . We prefer the Catalan number 
series z(x) over the modified Bessel function mainly because the inverse function of z(x) 
takes a simple form x = z + - . 

Motivated by the technique developed in (7J [Ml E7] for single Hurwitz numbers, let us 
define 

(7.17) = E (fcV^' 

fc=0 v y 



and 
(7.18) 

We then have 



, , ^/2fc + l\, _ 1 



(7.19) {„„(«) = -i + 1) 6.W = - i'] 6.W 

and 

(7.20) , n+l( t) = -± + 2) r^) = - 1 ) *,(*). 
The initial values are computed as follows: 

(7 - 21) ^ = 2( 1 "^J^^l(mJ^ 

1 / z(z 2 + 1) d 



2 V z 2 -\ dz z 2 - 1 4* 



and similarly 

(7.22) r/ (i) = - (/+ " 



At 

We note that £ n (t) and ?7n(i) are Laurent polynomials of degree 2n + 1 for every n > 0. 
Since they are defined as functions in x, we have the reciprocity property 

{ ' ' r, n {l/t) = r] n {t). 
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This follows from 



x 



-x. 



The annulus amplitude requires (g, n) = (0, 2) Gromov-Witten invariants. They can be 
calculated from the (g,n) = (0, 1) invariants using the Topological Recursion Relation [28J. 
The results are 



(7.24) (t w (u,)t m ») 



= f (mi!)^m 2 !) a (m 1+ m 2 +l) Ml = 2m U M2 - 2m 2 

I (m 1 !)2 1 (m 2 !) 2 ("1+1+2) ^1 = 2mi + 1, /J 2 = 2m 2 + 1. 



Theorem 7.7. T/ie annulus amplitude is given by 

(7.25) F^ 2 (z 1 ,z 2 ) = -log(l-z 1 z 2 ). 

Hence we have 

dt\ ■ dt 2 dx\ ■ dx 2 dt\ ■ dt 2 



(7.26) 



WUh,t 2 ) 



( tl -t 2 )2 ( Xl - X2 )2 (t 1+t2 )2' 



Proof. From (7.24) we calculate 



i °° 11 

^S^l,^) = X] ( T MiM r M 2 ( W ))o,2 l*lW -7TT+T -TT^+T 



MliA»2=0 



1 



2mA /2m 2 \ 1 



^ n (rrn + m 2 + 1) V mi y V ^2 y x? mi+1 x 2 ."^ 1 
mi,m2=0 i ^ 



1 



+ E 



mi ,m2=0 



(mi + m 2 + 2 



-(2mi + l)(2m 2 + 1] 



2mA /2m 2 \ 1 



m 1 )\m 2 J x \ mi+2 x 2 2 m2+2 ' 



Thus we have 



d d \ P i 



oo 

E 

mi,m,2=0 



2mA /2m 2 \ 1 



m l )\m 2 ) xl mi+1 x 2m2+1 
2mA /2m 2 \ 1 



2 E p™ 1+ i)p™ 2+ i)( 2 „:;)( 2 ™ 2 )^ 

mi ,m2=0 1 



1 



,2mi+2 2m 2 +2 



-2^0(^1)^0(3:2) - 2z'(xi)z'(x 2 ) 



-2 



x~l 



Z-2 



z\-\ zl-l 



X, 



-2 



ZlZ 2 (l + ZiZ 2 ) 

(z 2 -l)(z 2 2 -l) 



where £o(x) is calculated in (7.21), and from (4.6) we know 

~2 



z '( x ) = -jz 



dz 

dx z 2 — 1 



On the other hand, 



•''V; - ^l-M 



_ f Ziizf^rl) d^_ Z 2 {z\ + 1) _d 



z{ — \ dz\ 
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+ 



111 — Z\Z2 



ZlZ 2 (l + z l z 2 ) 



Therefore, 
d 



Xl ~dx[ + X2 ~d~x~ 2 ^ ( F Q> 2 ( Zl ' : " ! + l(>g! ' ~ l< '-■ 



d d 

Xl- h X 2 ~ r ~ 

ax i dx2 



\ ( °° 11 
J E (Vi(^H)o,2WW!^i^ +i 

00^/00 j 00 l\ n 

~ E n ( E ^ 2m+l E ^ 2m+l ] 
n=l \m=0 1 m=0 2 / 



0. 



Since the kernel of the Euler differential operator is the constants, and since actual compu- 
tation shows that the first few expansion terms of the Laurent series 

00 ^^00^/00 ^00 l\ n 

E < r /*»^»>0,2 HiW " E " E C -3^+T E C ™^+T 

/il,M2=0 X l X 2 n=l \m=0 X l m=0 X 2 / 



are 0, we complete the proof of (7.25). 



□ 



Using £ n (t) and rj n (t) of (7.17) and (7.18) and the classical topological recursion relation 
[28], we can systematically calculate the Laplace transform of stationary Gromov-Witten 
invariants. First let us consider (g, n) = (0,3). Since the sum of the descendant indices of 

(r w (w)T M (w)r w (w)) 0) , 

is even, we have 

(T2 mi {L0)T 2m2 {u)T2m 3 {u)) 03 ~ 

(7.27) 

(T2mi(w)r 2m2+ i(a;)r2 m 3 + i(w)) 3 



2 2 2 ' 



(m 2 + l)(m 3 + 1) 



m 2 (m 2 + l) 2 (m 3 + l) 2 
The Laplace transform is therefore 

(7.28) F^ 3 (t 1 ,t 2 ,t 3 ) = E ^i( w )^M T «M)o/i!f2!w!^ 



/.tl + l _M2+1 „./-»3+l 



E 

mi,m2,m3>0 



2mA f2m2\ (2mz\ 1 

"T-l / V m 2 



1 



1 



m 3 



)x\ 



2mi+l „2ro2+l 2m 3 +l 



+ E 

mi,m2,m3>0 



2mA [2m 2 + 1\ /2m. 3 + 1 



mi 



m 2 



m 3 



v - /2mi + A /2mA /2m 3 + A 1 1 

+ 2^ I m J [ m J y m J 2m 1+ 2 ■ 2m, 
mi,m2,rrt3>0 / \ / \ / y & 



+ E 

mi,m2,m3>0 



2mi + A /2m 2 + A / 2m 3 \ 1 



mi 



m 2 



m 3 



r 2 mi +l 


2m 2 +2 
x 2 


™2m 3 +2 


1 


1 


1 


2mi+2 


2m 2 +l 
x 2 


2m 3 +2 


1 


1 


1 


Jmi+2 


„,2m 2 +2 
x 2 


™2m 3 +l 
x 3 



&(*i)6(*2)6)(t3) + Co(*i)»»(*2)»»(*3) + r) (h)(o{t 2 )vo(t3) + Vo(h)vo(t2)Zo(t 3 



16 



(tl + l)(*2 + l)(*3 + l) 1 



^1*2*3 
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which is indeed a Laurent polynomial. Since it is an odd degree polynomial in £ n (i)'s, we 
have the reciprocity 

F^{l/t u l/t 2 ,l/h) = -i$(ti,*2,* 3 ). 
The n = 1 stationary invariants are concretely calculated in [66]. We have 

1 \ 2 / 1 , . , 1 /2d - 1 



(7.29) {T2d+2) ^-[d\) v.5!4^ (2(i " 1) + 242 V 2 

(^-^.i = (£) e( 2 ^ x ) e np^p-' 

£ — 1 fcj ^0 i — 1 
felH hfcf=9 

We thus obtain 

(7-30) = 1 f f") - 1)^T = M " 6 



b(*i)) 



1 / « _ 7 1 



384 V" 1 7h + *i if 



To calculate the g = 2 case we need to do the following. 



„2d+3 



(t 2 -!) 2 d 



„2d+l 



8i di 
1 (t 2 - l) 3 



5j^(2^i(*i) " &(*!)) + ^ ( 2 &(tl) " 36 (*i) + &(ti)) 
525if - U70t\° + 1107*? + 527if + 1107^ - 1470t 2 + 525). 



2 19 • 3 2 • 5 t 9 

Proposition 7.8. Fj^ti) is a Laurent polynomial of degree 6g — 3 with the reciprocity 
Proof. First we calculate the binomial coefficient 



(2d - l)(2d - 2) • • • (2d -e) = jj [2 l d l - v 2 J d l ~ l + ■■■ + {-1)H\ 

as a polynomial in d, and then replace each d l with £i(ti). The result is a linear combination 
of £o(*i)> • • • iii{t\)- Let S^(ti) denote the resulting Laurent polynomial of degree 2£ + 1. 
Then we have an expression 

fciH — hke=g 



40 



O. DUMITRESCU, M. MULASE, B. SAFNUK, AND A. SORKIN 



it" 



8t 



dtj 



2g-2 



e s«<«o e n 



ki>0 
felH h^=9 



L (2fc + 1)! 4fr 



which is a Laurent polynomial of degree 2{2g — 2) + 2g + l = 6g — 3. The reciprocity property 
follows from (7.23) and the x\ expression of (7.32), where x\ changes to —x\. In particular, 
the even order differentiation in x\ is not affected by this change. □ 

The Eynard-Orantin recursion is for the differential forms Wj n (ti, . . . , t n ). We need the 
recursion kernel. From ( |7.13 ) and (7.26), we compute 



(7.33) K ¥ (Ml 



2 Wg(-t)-Wji(t) 



t + ti + i - h 



(t 2 



1 



io s I (trip 



l 



l l 

+ 



16 + * — *i 

We note the reciprocity property of the kernel 

(7.34) K pl (l/t,l/t 1 ) = -K ¥ \t,t 1 ). 



-8tdt 
1 



dti 

2 



log v (*+i) 2 



The topological recursion (2.3) becomes 

i^ pl (Mi) 



(7.35) W^ n (t u t 2 ,...,t n ) 



1 

2^i 



Wg_ 1>n+1 (t, ~t,t 2 , ...,t n ) 



No (0, 1) terns 

+ E <i/i + i(*>*/)<Vi + iH>^) 

91+92=9 
7Uj={2,3,...,n} 



where the residue calculation is taken along the integration contour 7 (see Figure 4.1) 
consisting of two concentric circles of radius e and 1/e for a small e centered around t = 0, 
with the inner circle positively oriented and the outer circle negatively oriented. Since 
there is a log singularity in the complex i-plane, we cannot use the residue calculus method 
to evaluate the integral at t = t\ and t = —t\. Thus the residue calculation of (7.35) is 
performed around the neighborhood of t = and t = 00. 

So let us provide two expansion formula s for the kernel K p (t, t\), assuming that t\ € C* 
is away from the log singularity of Figure 



7.1 



The transcendental factor of K (Ml) nas 



an expansion 
(7.36) 



1 1 



44 



428 



f]n(r ((*-i) 2 



I I t J f* 4. 

t 2 3 45 945 14175 



t° + 



10196 

467775' 



:t° + ■■■ 



around t = 0. The denominator of the coefficient of t 2k 2 is given by 

n q l^tJ = 3 l^ -5LU -tlij ... , 

5=3, prime 
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which is the same as ^(L/.) of |38^ Lemma 1.5.2]. The expansion of 
given by 



+ 



at t = is 



+ 



t + ti t-h 



-2t 



1 1 



t 2 

71 



- 2 E 



f 

=0 l l 



2n+2 ' 



From the expression (7.33) and the above consideration, we know that around t = 0, 
K^ 1 (t, t\) starts from t" 1 , and that the coefficient of t 2n ~ l is a Laurent polynomial in if 



starting from 
(7.37) K v \tM) 



(2n+2) 



up to t-y with rational coefficients. More concretely, we have 



1 1 

32 if 



+ t 



1 1 

32 tf 



7 1 

96 if 



+ t d 



+ t a 



1 1 

32? 



96 i? 1440 30240 tf 



+ t' 



1 1 

32 tf 



7 1 71 1 

+ 

1 UK} b i 
7 1 71 1 191 

96 tf + 



1 1 

32^ 

191 ] 



7 1 71 1 

96 tf + 1440? 



1 



23 1 



+ t M 



7 1 



96t}° 



+ 



1 1 

.32 tp 

Similarly, around t = 00 we have 

1 



71 1 



1440 tf 30240 tf 28350 tf 
191 1 23 1 233 1 



1440 tf 



30240 tf 



28350 tf 935550 if 



+ 



dt 



(7.38) K r (t,h) 



t" hi 

32 



1 2 7 

— i? H 

32 1 96 



32 



96 



71 
1440 



+ 



1 

+ ? 



32 



32 1 96 1 



71 2 
*1 + 



1440 



191 \ 
30240 J 



+ 



t 7 



32 1 96 1 



i? + — i 
1 96 

71 



6 



1440 



71 
1440 

191 
30240 



191 



30240 
23 



23 



28350 



28350 



233 \ 



935550 



dt 



Theorem 7.9. The Eynard-Orantin differential form Wg n {t\, . . . , t n ) is a Laurent polyno- 
mial in tf , t|, . . . , t n 0/ degree 2(3g — 3 + n) in the stable range 2g — 2 + n > 0. /t satisfies 
the reciprocity property 

(7.39) Wf»(V*l, • • " > l/tn) = ("l) n <n(tl, • • • , tn) 

as a meromorphic symmetric n-form. The highest degree terms form a homogeneous poly- 
nomial of degree 2(3g — 3 + n), which is given by 

2k 1 



(7.40) WUh,. 



■ j tn) 



22g-2+2r 



^2 ( T ki--- T kJg,nY]_ 



,k n >0 



i=l 



(2h + 1) 



dt; 



Indeed it is the same as the generating function of the ip-class intersection numbers (5.4)- 



Proof. The statement is proved by induction on 2g — 2 + n using the recursion (7.35). 
The initial cases (g, n) = (1,1) and (g,n) = (0,3) are easily verified from the concrete 
calculations below. Since we are expanding j^- + j^j- around t = and t = 00, it is 
obvious that the recursion produces a Laurent polynomial in if , t|, . . . , i^ as the result. 
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The expression of (7.38) tells us that the residue calculation at infinity increases the 
degree by 4. This is because the leading term of the coefficient of t _ ( 2n + 1 ) is if +4 , and the 
residue calculation picks up the term t 2n . By the induction hypothesis, the right-hand side 
of ( 7.35 ) without the kernel term has homogenous degree 2(3g — 3 + n) — 4. The reciprocity 
property also follows by induction using (7.34). 

The leading terms of Wg n (t\, . . . , t n ) satisfy a topological recursion themselves. We can 
extract the terms in the kernel that produce the leading terms of the differential forms from 
Q7.36H or ( |7.38D . The result is 



(7.41) K WK (t,h) 



-t 3 y 

32 ^ 



k=0 



if 1 
l^di 



t - U 



+ 



t + ti 32 



dt 



which is identical to [TO], Theorem 7.4], and also to (5.10). Since the topological recursion 
uniquely determines all the differential forms from the initial condition, and again since the 
(g,n) = (0,3) and (1,1) cases satisfy ( |7.40 ), by induction we obtain (7.40) for all stable 
values of (g, n). □ 



The (g, n 
(7.42) WfUh) 



(1, 1) Eynard-Orantin differential form is computed using (2.4) 

dx ■ dx\ 



1 






S K 


" 2vri 




if 




64 I 


2lTiJ 




dt ■ dt 



dt dt 



384 1\ 



1 1 

128 if 



dU 



This is in agreement of Wfi(ti) = dF^^i) and (7.30). From (7.35) we have 



(7.43) Wl 3 (h,t 2 ,t 3 ) 
1 

~ 2tH 



K v \t,h) [<2(i,*2)<2(-i,i3) + <2(i,i3X2(-*,i2) 



1 

16 



1 + 



j-2 j-2 j-2 

l 1 L 2 (-3 



dt\dt2dU 



It is also in agreement with (7.28). 

Norbury and Scott conjecture the following 

Conjecture 7.10 (Norbury-Scott Conjecture 



For (g, n) in the stable range we have 



(7.44) 



Wq tn (h, ■ ■ ■ ! t n ) — dl ■ ■ ■ d n Fy (tl, . . . ,t n ). 



The conjecture is verified for g = and g = 1 cases in [62] . We recall that the Eyanrd- 
Orantin recursion for simple Hurwitz numbers is essentially the Laplace transform of the 
cut-and-join equation [23]. For the case of the counting problem of clean Belyi morphisms 



the recursion is the Laplace transform of the edge-contraction operation of Theorem 3.3 



Question 7.11. What is the equation among the stationary Gromov-Witten invariants of 
P 1 whose Laplace transform is the Eynard-Orantin recursion (7.35)? 
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Appendix A. Calculation of the Laplace transform 



In this appendix we give the proof of Theorem 4.3 



Proposition A.l. Let us use the Xj-variables defined by Xj = e Wj , and write 

W® n (h, ...,t n ) = w g , n (x\, ... ,x n ) dx\ ■ ■ ■ dx n . 
Then the Laplace transform of the recursion formula \3. 15 ) is the following differential 



recursion: 



(A.l) -x 1 Wg 7n (x 1 ,...,x n ) 



E 

J'=2 



d_ 
Ox-, 



(Wg,n-l(%2, ■■ ■ ,X n ) ~ Wg,n-l(xi, %2, ■ ■ ■ , Xj , . . . , X n )) 



+ W g - ljn+1 (x 1 ,X 1 ,X 2 ,...,X n )+ ^ w g 1 ,\I\+l(xi,Xl)w g2j \j\ +1 (x 1 ,Xj). 



91+92=9 
IUj={2,...,n} 



Proof. The operation we wish to do is to apply 

n 

(-1)" £ ^...^JJ-s+r 



A t lv,Atn>0 



i=l X i 



to each side of (3.15). Then by (4.13), the left-hand side becomes w g ^ n (xi, . . . ,x n ). 

The second line of (3.15) is straightforward. Let us just consider the first term, since the 
computation of the second term is the same. 

- 1 



i=l a+/3=tM 1 -2 



(_l)«+i ^ a/3/x 2 ---/in^-i,n+i(a,/3,/i2...,^n)-^ TT • ^rn^+T 

1 At2,-^n>0a,/3>0 X l X l i=2 X i 



x\ 



W g -l, n+ l{xi,Xl,X 2 , ■ ■ ■ ,X n ). 



Thus the second line of (3.15) produces 

-—{ w g-i,n+i(xi,xi,x 2 ,...,x n )+ ^ 

Xl V 91+92 = 



w Sli |/| + i(a;i, xi)w 92i \ j\+i(xi, xj) 



IUj={2,...,n} 



To calculate the operation on the first line of (3.15), let us fix j > 1 and set v 
Hi + /J,j — 2 > 0. Then 

(A.2) (-l) n J2 M2---Mn(w + Mi-2) 

[ll,...,fj, n >0 

n 1 

X D g ^i{m +flj - 2,^2, ■ ■ ■ ,Ttj, ■ ■ ■ ,/Un) 

1=1 x j 

oo 

= -J] Yl {-l) n ~ l VH2---'fc---Hn 

v=0 H2,—4lj,—,Hn>0 
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1 1 U+1 

X £> fl ,„_l(l/,/Z 2 ,... > ^V--,A*n)-£+r ]J Yl ^ X l 3 

Xi ■ i ^ ■ X- 



-2 l 



1 i^hj X 



M+ 1 



Assuming < \xj\, we calculate 

v+l 1 -t n v+l 



(A.3) Wi j 



-2 1 



1 d 



M+ 1 



xi y s 1 9 



( 



1 



i/+2 1 



1 _ £l 1 _ £1 



1 .9 / 1 



<9xj \ 1 



£1 

Xj 



+ x\ 



1 1 



5 



dxj \xj - xi 



v+l ■ 



We then substitute (A.3) in (A. 2) and obtain 



1 d ( 1 

(A. 4) (A. 2) = w g , n -i(xi,x 2 , ...,£j,.. .,x n )-2 ' — " 



£1 

Xj 



x\ dxj 



x{ dxj \ 1 

\Xj X\ J 



1 d ( 1 



x \ dx j \^ x j x \ 
This completes the proof. 



\Wg,n— 1 i^2i • • • j ICj i ■ • • > Wg,n— 1 (^1 1 X2 , • • • , 37 j , • ■ • , %)) I ■ 



□ 



Proof of Theorem J±.3. When the curve is split into two pieces, the second term of the third 
line of (A.l) contains contributions from unstable geometries (g,n) = (0, 1) and (0,2). We 
first separate them out. For gi = and / = 0, or g<i = and J = 0, we have a contribution 
of 

2wo t i(xi)w g>n (xi, X2, ■ . -,x n ). 
Similarly, for g\ = and / = {j}, or g 2 = and J = {j}, we have 

n 

(x\, . . . 

3=2 

bmce and W ° 2 are defined on the spectral curve, it is time for us to switch to the 
preferred coordinate t of (4.1) now. We thus introduce 



(A. 5) W gn (ti, ...,t n ) = w gn (ti, ...,t n ) dti ■■■dt n = w g>n (xi, ...,x n ) dxi ■ ■ ■ dx n . 
Since WQ t i(x) = —z(x), we have 

wo,i(x) 

WQ,2(xi,X 2 ) 



t+1 
t-1 

1 (t\ - l)^ if 2 - I) 2 



(h + h) 2 8*i 8t 2 

n U2 _ 1 \2 

W g>n (x 1 ,...,X n ) = (-l)Xjil,..,t„)Il ' , ■ 



i=l 



Thus (|A.1|) is equivalent to 
2 



t\-\ tl _ 1 J W 9,n( t ^--^ t n) 
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" ' (*j-l) 2 (*?-l) 2 8tj d . r , . 



+ 



a / ( tf- 1)^-1) gti (tf- i) 2 D 

dtj { 4(tj - t 2 ) {t\ - l) 2 8t 



UVn-l(*2,---,tn) 



+ 



8ti 



stable 



\ 



l,n+l(*l5*l,*2, • ■ ■ ,t n ) + £ 

IL 

1 (*? I) 2 .. 



71' 



,|J| + l(*l>*/)w£,|./| + l (*!>*./) 



91+32=9 
/Uj={2,...,n} 



= E 

i=2 



i=2 

+2 i ^2 



(tl+tj) 2 8 *1 9 



^■(tf-1) 2 , 1 (tf-1) 2 



+ 



W o,ri-l(^l' ■ ■ ■ ,tj, ■ ■ ■ , tn) 



ti s / (t 2 -i) 3 „ 



+ 



q-i dtj utjitj-m 



2 _ 1 2 A "'</•" M'-'- ■ • ■ >*ra) 



+ 



(*? ~ I) 2 

8ti 



stable 



\ 



,D 



W, 



l t n + l(t 1 ,t 1 ,t 2 ,...,t n )+ £ w f 1 ,\I\+l( t ^ t l) w f 2 ,\J\+l( t U t j) 



91+92=9 
IUj={2,...,n} 



J 



Since 



we obtain 

n 

(A.6) w° n (t 1 ,...,t n ) = -Y, 

3=2 



t\ + i _ h + 1 \ _ _ 4ti 
— l h-i) ~ tj-V 



^1 16^2— ^ «W(*2,...,tn) 



+ 



16tf (t 2 - ^ 2 9 



^Vn-l^l) • • • , ij, • • • , in) 



(*? ~ I) 3 

32t 2 



stable 



l,n+l(*l»*l)*2)---,*n)+ E ^,|/|+l(*l'^)^,|J|+l(*l'^) 



V 



Now let us compute the integral 

1 1 



(A.7) W? n (t 1 ,...,t n ) = 



91+92=9 
IUj={2,...,n} 



I 1 \ (t 2 - l) 3 1 , 



64 2ni L \t + ti t-hj t 2 dt 



n , 

£ I W,g(f, *2, • • ■ , tj, . . . , t n ) + W£ 2 (-t, t^Wg^it, t 2 , . . . , tj, . . . , t n ) 

j=2 V 



46 



O. DUMITRESCU, M. MULASE, B. SAFNUK, AND A. SORKIN 



stable 



+ w, 



D 

g-l,n+ 



1 (t,—t,t2,---,t n )+ 



w. 



D 

91,1*1+1 



31+92=9 
IUj={2,3,...,n} 



Recall that For 2g— 2+n > 0, Wg n {t\, . . . , t n ) is a Laurent polynomial in q, , 
third line of (A.7) is immediately calculated because the integration contour 7 of Figure 4.1 



(M/^iji+^-Mj) 



, t 2 n . Thus the 



encloses ±t\ and contributes residues with the negative sign. The result is exactly the last 
line of (A.6). Similarly, since 

W D 2 (i, tj)W g)n - X {-t, t 2 , . . . , tj, . . . , t n ) + W D 2 (-i, ij)W fl , n _i(t, i 2 , • . • , tj, . . . , t n ) 
1 1 



+ 



. 3 , (t-tj) 
the residues at ±t± contributes 

4-2 i\3/+2 



W 



D 



n _l(i, t2, ■ ■ ■ ,tj , . . . , t n ) dt dt dt 2 ■ ■ ■ dtj ■ ■ ■ dt n , 



,tj 1 ... ) t n ). 



This is the same as the second line of the right-hand side of (A.6). 

Within the contour 7, there are second order poles at ±tj for each j > 2 that come from 
WiP 2 (±£, Note that acts as the Cauchy differentiation kernel. We calculate 



1 1 

64 2tH 



+ 



t + ti t-ti 



(t 2 



t 2 



i\3 n 

— E 



i=2 



+ w 0,2(- t > t j) w g,n-l(t,t2 



(t] ~ 1) 





This gives the first line of the right-hand side of (A.6). We have thus completes the proof 
of Theorem 14.31 □ 
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